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PREFACE 
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NOTATION 


a negative of angular momentum per unit mass of the rotating central body 

A{u) linear function of u defined by Equation 48 

A*{u *) linear function of u* defined by Equation 19 

B(u) cubic function of u defined by Equation 45 

B*(u*) cubic function of u* defined by Equation 20 

C*(u*) fifth degree polynomial in u* defined by Equation 21 

D(u) second degree polynomial in u defined by Equation 46 

D*(u*) second degree polynomial in u* defined by Equation 27 

e constant of the motion related to the energy of the particle, given by Equation 38 for 
equatorial orbits and by Equation 89 in general 

e* constant of the motion related to the energy of the particle, given by Equation 24 for 
equatorial orbits and by Equation 1 22 in general 

E{u) linear function of u defined by Equation 47 

E*(u*) linear function of u* defined by Equation 28 

f(u) function defined by Equation 165 

g the metric tensor 

h constant of the motion related to the angular momentum of the 
tion 37 for equatorial orbits and by Equation 88 in general 

h* constant of the motion related to the angular momentum of the 
tion 23 for equatorial orbits and Equation 1 2 1 in general 

H quantity defined by Equation 156 

m mass of the rotating central body 

r radial coordinate in S form defined by Equation 1 5 

r* radial coordinate in E form defined by Equation 8 

s proper time 


particle, given by Equa- 
particle, given by Equa- 
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t 

t* 

u 

u* 

U B1,B2,B3 

f 

U B1 ! B2 


U Dl,D2 

U E 

Uf 

U s 

u \, II, III 
u Lli 

V 

W 

x,y,z 

P 
A 
A* 
6 , $ 
d*, <p* 


a time coordinate in S form defined by Equation 1 5 
a time coordinate in E form defined by Equation 8 
inverse radius, 1 /r 
inverse radius, 1 /r* 

roots of the cubic B(u ) = 0, which are also denoted by Wj u m in the figures 
roots of the quadratic B'(u) = 0, which are also denoted by u' x n in the figures 

root of B’\u) = 0 (also denoted by u s ), which is the point of symmetry of the cubic 

B(u) 

repeated root of f(u) = 0 

roots of the quadratic D(u) = 0, given by Equation 197 
root of the linear equation E(u) = 0, given by Equation 206 
root of /(u) = 0 
see 

See U B1,B2,B3 
see u BX B2 

quantity defined by Equation 141 

a measure of the energy of the particle, given by Equation 153 

spacelike coordinates in M form defined by Equation 1 

quantity equal to 2m 

quantity defined by Equation 14 

quantity defined by Equation 69 

angular coordinates in S form defined by Equation 1 5 

angular coordinates in E form defined by Equation 8 
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A QUALITATIVE STUDY OF THE COMPLETE SET OF SOLUTIONS 
OF THE DIFFERENTIAL EQUATION OF MOTION OF A TEST PARTICLE 
IN THE EQUATORIAL PLANE OF THE KERR GRAVITATIONAL FIELD 


by 

Harry E. Montgomery 
Goddard Space Flight Center 
and 

F. K. Chan* 

Catholic University of America 


INTRODUCTION 

In 1963, Kerr (Reference 1) presented an exact solution of Einstein’s empty-space field equations 
that might be the exterior field of a rotating body. He presented his metric in two sets of coordinates 
which were later labeled as the E and the M coordinates. He also stated that the two parameters in his 
solution are the Schwarzschild mass and the specific angular momentum of the central body. 

In 1964, Boyer and Price (Reference 2) used Kerr’s metric in E coordinates and obtained the equa- 
tion of motion of a particle in the equatorial plane and then studied the precession of the pericenter of 
the orbit. By comparing this result with the motion of the pericenter as determined by Lense and 
Thirring (Reference 3), who had obtained in 1918 the perturbations of the elements of planetary orbits 
in an approximate gravitational field of a uniformly rotating sphere of constant density, Boyer and 
Price identified the two parameters m and -a in the Kerr metric as the mass of the sphere and the spe- 
cific angular momentum of the sphere, thus confirming Kerr’s assertion. 

In 1 966, Carter (Reference 4) studied the motion of a particle along the symmetry axis of the 
Kerr gravitational field using E coordinates. 

In 1967, Boyer and Lindquist (Reference 5) discussed some coordinate systems useful in the study 
of the Kerr metric which they called the Minkowski-like M, Eddington-like E, and Schwarzschild-like S 
coordinates. Certain properties of the Kerr metric, including the behavior of geodesics lying in the 
equatorial plane, were examined in some detail. 

In 1968, DeFelice (Reference 6) investigated the geodesic motion in the equatorial plane of the 
gravitational field of a rotating source using E coordinates. The differential equations of the geodesics 
were studied to obtain the shape of the loci of turning points of test particles without introducing 
numerical values of the parameters involved. Figures were given that show the behavior of geodesics 
for different ranges of these parameters. 

*F. K. Chan is presently with Computer Sciences Corp., Silver Spring, Md. 20910. 
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In 1968, Carter (Reference 7) demonstrated that the Hamilton-Jacobi equation for a test particle 
using E coordinates can be separated. 

The present work deals with the motion of a particle in a Kerr gravitational field. It begins with 
a brief review of the Kerr metric in the three coordinate systems and also a derivation of the equations 
of motion of a test particle. These equations are derived in the E and the S coordinate systems. 

Next, a step-by-step procedure is presented for the qualitative analysis of the complete set of 
orbits of a test particle in the equatorial plane, using S coordinates. The results for each value of a are 
represented in a master diagram that gives the complete topological genera of the orbits for various 
ranges of the two parameters: the angular momentum per unit mass of the test particle h and the 
“energy” per unit mass of the test particle e. In addition, the physical orbits corresponding to the 
master diagram are given. 

A discussion is then given of the determination of the direction of orbital motion, the scaling of 
solutions for j3 ¥= 2, conversion of results of the present S coordinates to E coordinates, and, finally, 

the p seudo singularity at u = u D . 

The method of solution used in the analysis relies heavily upon Chan’s qualitative solution of the 
motion, of a particle in the Schwarzschild gravitational field. Because this work is unpublished, it is 
presented in Appendix A. The algebraic details that were omitted in the text are included as Appen- 
dix B. 

EQUATIONS OF MOTION OF A TEST PARTICLE IN THE KERR GRAVITATIONAL FIELD 


Various Forms of the Kerr Metric 


For a rotating body, the exterior gravitational field is described by the metric originally given by 
Kerr (Reference 1) in the form (henceforth referred to as the Minkowski-like M form) 


2mr* 2 k 2 

ds 2 = -dx 2 - dy 2 - dz 2 + dt* 2 , 

r * 4 + a 2 z 2 

where k is given by 

(. r * 2 + cc)r*k = r* 2 (x dx + y dy) + ar*(x dy - y dx) + (r* 2 + a 2 )(z dz + r* dt*) . 


( 1 ) 

( 2 ) 


The coordinate r* is defined by 
where 


r* 4 - {R 2 - a 2 )r* 2 - a 2 z 2 = 0 , 
R 2 =x 2 + y 2 + z 2 . 


(3) 

(4) 


In these equations m can be interpreted as the Schwarzschild mass and -a the angular momentum per 
unit mass about the z axis. However, in this form of the line element, it is difficult to understand truly 
the physical nature of the “radial” coordinate r* because of the presence of too many off-diagonal 
terms in the metric tensor g, iy . Moreover, setting the parameter a equal to zero does not simply yield 
the well-known Schwarzschild metric. In his work, Kerr also introduced the following transformation, 

x = {r* 2 +a 2 ) 112 sin0* cos - tan -1 ~~j, (5) 
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y = (r* 2 + a 2 ) 1 / 2 sin d* sin (<p* - tan 1 — 

V r*i 


( 6 ) 


z = r* cos 6 * , 

to obtain the line element in the form (henceforth referred to as the Eddington-like E form) 
ds 2 = -dr* 2 - 2 a sin 2 6* dr * dip* - (r* 2 + a 2 ) sin 2 6* dip* 2 - P* dd* 2 + dt* 2 

- (dr* + a sin 2 6 * dip* + dt*) 2 , 


where P* is defined by 


p* = r* 2 + a 2 cos 2 6 * . 


(7) 

( 8 ) 
(9) 


Even now with the elimination of some off-diagonal terms in g , it is still difficult to interpret the 




radial coordinate r* or to obtain the Schwarzschild metric directly by setting a - 0. 

However, Boyer and Lindquist (Reference 5) later introduced the following transformation: 


r-r*. 

(10) 

6 = 6* , 

(11) 

a dr* 

dip = dip * + , 

A 

(12) 

Br* dr* 

dt = dt* - — - — — , 

(13) 

A = r 2 - Br + a 2 , 

(14) 


where 

which then reduces Equation 8 to the form (henceforth referred to as the Schwarzschild-like S form) 

\2 


( dr 2 \ Br 

— + dd 2 ) - (r 2 + a 2 ) sin 2 6 dip 2 + dt 2 - — (a sin 2 6 dip + dt) 2 


where 


P = r 2 + a 2 cos 2 6 . 


(15) 

( 16 ) 


The outstanding feature of this equation is that it has only one off-diagonal component g 0f and reduces 
nicely to the standard Schwarzschild metric when a = 0. 

The transformation given by Equations 10 to 14 reveals that there is a mathematical pseudo- 
singularity introduced when A = 0. Even though in this respect the E coordinates appear to be better 
for analysis of the Kerr metric, it turns out that the S coordinates are preferable in terms of physical 
significance in that they determine the location of horizons or stationary null hypersurfaces when 
A = 0. (For a more detailed discussion, see Boyer and Lindquist (Reference 5).) 
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Equations of Motion for a Test Particle in the Equatorial Plane 


Boyer and Price (Reference 2) used the E form of the Kerr metric in their study on the interpre- 
tation of the parameter -a as the angular momentum per unit mass of the central spinning body. For 
test particles with orbits in the equatorial plane they obtained the following equation of motion: 

/ du *\ 2 du* 

A *(u*) — +2 aB*(u*) C*(u*) = 0, (17) 

\d<t>*} y 'd<p* 

where 

u*= — , (18) 

A*(u*) =-(3(h* + e*a) 2 u* + h * 2 + (1 - e* 2 )a 2 , (19) 

B*(n*) = (3(h* + e*a) 2 u*^ - [h* 2 + (1 - e* 2 )a 2 ]u* 2 + (3u* - (1 - e * 2 ) , (20) 

C*(u*) = (3a 2 (h* + e*a) 2 u* 5 - ||3 2 (/z* + e*a) 2 + a 2 [h* 2 + (1 - e* 2 )a 2 ] } «* 4 

+ 2 (3[h*(h* + e*a ) + a 2 ]u* 2 - [h* 2 + (3 2 + 2(1 - e* 2 )a 2 ]u* 2 

+ j3(2 - e* 2 )u* - (1 - e* 2 ) , (21) 


(3 = 2m , 


h* — a ( 1 + — ) r* + ( r* 2 4 - a 2 + 


Pa : 


Ba . 

6* + — t* 


7 ^ 


P Pa . 


I* 5 }* 


*V" 

It can be shown that 


0* + 1 T t* . 


where 


C* =B*D* 
A*D* + a 2 B* =E* 2 , 
D*(u*) = a 2 u * 2 - 1 3u* + 1 
E*(u*) = ~(3(h* + e*a)u + h* . 


( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


Next, the equation of motion in the S form will be derived from first principles. (It may also be 
obtained from Equation 17, which is for the E coordinates, and Equations 10 and 12, which are the 
equations of transformation from (r*, </>*) to (r, <p). A method of finding equations of motion corre- 
sponding to a general metric is given by Synge (Reference 8). If the line element is denoted as 


ds 2 = S M! , dx^ dx v 


ix,v=\, 2, 3, 4 


and the Lagrangian F(x, x) is defined by 


F(x, x) = x v , 


(29) 


(30) 
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then the equations of motion are given by 



— — - — =0 for ti= 1,2, 3,4, 

* 3i M 3x M 

( 51 ) 

where 

• M — 

X ds 

(32) 

Now, 

for the case of equatorial orbits we have 6 = ir/2, so that from Equation 1 5 the Lagrangian 

becomes 

y 2i2 « 

F = (r 2 + fl 2 ) 0 2 + i 2 - - {a 0 + t) 2 . 

A r 

(33) 

Hence, the 

following three equations corresponding to r, 0, and t result: 



rf 3F 3F _ o 
ds JF 3 r " ’ 

(34) 


d 3F 3F _ 
30 30 

(35) 


d 3F 3F _ o 

ds Pi 3 t 

(36) 


These three equations and the line-element equation (Equation 1 5) constitute a set of equations of 
which only three are independent. It is simplest to use Equations 15, 35, and 36 to obtain the equa- 
tion of motion in the equatorial plane. 


Because F, which is given by Equation 33, is not a function of 0, Equation 35 implies that 
3F/30 is equal to a constant -2h; therefore, 

3 F ~ ~ . 2)3 a 

— r = - 2{r 2 + a 2 ) 0 - (a 0 + t ) = - 2h . 

30 r 

Similarly, because F is not a function of t, PF/di is equal to a constant 2e; thus 

3 F .2 p ■ . 

— r =2 t - — (a 0 + f ) = 2e . 

3 1 r 

From Equations 37 and 38 for 0 and i, it follows that 

0 = -7 [h(r- P) - Pae] , 
rA 


t- — [e(r 3 + a 2 r + pa 2 ) + pah] . 
rA 


lfu= Mr, then 

Equations 39 and 40 then become 


•2 _ 0 2 MA 2 

m 4 \<i0/ 


(37) 


(38) 


0 = — [~|3(/z + ea)u + h ] , 
A 


(40) 

(41) 

(42) 
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u 2 A 


e(l + a 2 u 2 + (3a 2 u 2 ) + fiahu 2 


(43) 


Substitution of these equations into the line element given by Equation 1 5 yields 

/du\ 2 D 2 (a) 

— = B ^y 


\d<p ) 


E 2 (u) 


where 


B(u ) = j 3(h + ea) 2 u 2 - [h 2 + (1 - e 2 )a 2 ]u 2 + (3u - (1 - e 2 ) , 
D(u) = a 2 u 2 - (3u + 1 , 

E(u ) = -j3(/z + ea)w + h . 

In addition, A(u) is defined by 

A(u) = -j 3(h + ea) 2 u + h 2 + (1 - e 2 )a 2 . 

It follows from Equation 26 that the following identity holds: 

AD + a 2 B=E 2 . 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


Equation 44 is studied extensively in the next section and was chosen in preference to Equation 
17 because it is simpler and its metric reduces to the standard Schwarzschild metric when a = 0. 

With the S coordinates the equation of equatorial orbital motion has been obtained in a simpler 
form as given by Equation 44. The E form obtained by Boyer and Price as given by Equation 1 7 may 
also be reduced to a simpler form, as follows: Completing the square in Equation 1 7 results in 

aB*\ 2 


du* 


B* 


(A*D* + a 2 B*) = 0 . 


\d<p* A*/ A* 2 

By using the identity in Equation 26, Equation 50 takes the form 


/du* aB*\ 
\d<p* A*/ 


oB*\ 2 B*E* 2 

A* 2 


0 . 


Letting 


du + 

d<p + 


du* aB* 
+ 

d(j)* A* 


then 


u = u 
2 


'du + ' 

\d<f, 


B*E* 2 

A* 2 


(50) 

(51) 

(52) 

(53) 

(54) 


When a - 0, this equation also reduces to the equation of motion in Appendix A except for a factor of 
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h* 2 . Because u + = u*. Equation 52 yields 

A* d<p* dr* 

d(j> + = , (55) 

A * dr* - aB*r* 2 d(p* 

Of 

A * d(j) + dr* 

d(j>* = . (56) 

A* dr* + aB*r* 2 d(j) + 

Hence, Equations 53 and 55 give the transformation from (r*, <p *) to (r + , 0 + ). However, one 
drawback of this transformation is that Equation 55 involves the constants of motion h* and e* of the 
test particle because of the presence of B* in the equation. This is to be contrasted with the trans- 
formation from (r*, <p*) to (r, 0) given by Equations 10 and 12 for which no such weakness occurs. 


Equations of Motion of a Test Particle Not Restricted to the Equatorial Plane 


The differential equations for general orbits using the E form of the Kerr metric have been derived 
in this study. Because the method of derivation is the same as that used to obtain Equations 102 to 
105 in the S form, the resulting equations are presented, without derivation, for future reference. 

Using the following definitions. 


<• ^ — y^C' 


x| = 6* 


3 * ’ 

X* = t* , 


x* = r* 


x* = 6* , 


and letting the dot denote d/ds, the differential equations of motion are 


: v* 
X 5 


v* = V* 
x 2 6 


-l 


V * — 

x 3 


P* A* sin 2 xj 


(fiae*x* sin 2 x| + fih*x* - h*P* + aP*x% sin 2 x|) 


(57) 

(58) 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 

(65) 


— 

X 4 


jp* A* 

h* 

°i 


[ pah*x* + e*P*(x^ 2 + a 2 ) + e*pa 2 x* sin 2 x| + 0x*P*x£ ] , 


( 66 ) 

( 67 ) 


l 21 
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( 68 ) 


where 


a *i b *- a n b * 


n* a* 
U 12 a 21 


A* =x* 2 - (lx* + a 2 , 

2*, = 2x* + 2 ax*® 2 sin 2 x* + 2 a <j>* sin 2 x* + 2(x* 2 + a 2 ) 0*<i> 2 sin 2 xj 


2 P* (t* - e*)T ' 


2 t*T 2 + 


a* 2 = 2P* X * , 


a 21 =2 + 2a<& 2 sin 2 x| + 2T 2 


b f — 2&x|<I> 1 sin 2 x| - 4#x|x| 0* sinx| cosx| - 2x^x| 0* 2 sin 2 x| 

- 2(x* 2 +a 2 )x* 0* 2 sinx| cos x* - 2(x* 2 + a 2 )<$ 1 0* sin 2 x|-P*x* 2 +2 t*T x 

P* (/* - <?*) 2P* (/* - e*)P P*x* (/* - e*) 2 


b* — 2a$> x sin 2 x| - 4 a 0*x| sinx| cosx| - 2T X + 2xf 0* 2 sin 2 x| + 2x*x| 2 

PQ* 2 2 (3x* 2 Q 2 


P* A* 


- [ pah*x f + e*P*{x^ 2 + a 2 ) + e*pa 2 xf sin 2 x| + ^x^P*x^] , 


P* A* sin 2 xf 
P* . 

}* = _ e *) > 


(Pae*x* sin 2 x? + j 3h*x* - h*P* + aP*x? sin 2 x.?) , 


P* = 2xfxf - 2a 2 x5 


cos XT sm x: 


P* A* 


[/ pah*x% + e* P*(x^ 2 + a 2 ) + 2 e*P*x^x^ + e*j3a 2 x| sin 2 x| 


+ 2e*pa 2 x*xg sinx| cosx| +j 3P*x| z + 0x^P*x|] 


* 4 - 2 


p*A* 


[A* P* + P(2x^x| - £x|)] 
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( 80 ) 


Px* 

A* 


$ 


l 


$ 


2 


| -)3ae*x| sin 2 x| - 2j3 ae*x*x| sinx| cosx| 

P* A* sin 2 xj 

- )3/z*x| + h* P* - a P*x | sin 2 x| - 2aP*x|jc| sin xj cosx| 

- 0* [(2x*x| - j8x|)P* sin 2 x^ + P* A* sin 2 x| + 2P*A*x| sin x| cos x|] | , 

_1_ 

A* 


(81) 

(82) 


Equations 63 to 68 with appropriate initial conditions x*(s 0 ) = x? can be solved by numerical 
integration. Only five of the values of xf may be chosen independently and the remaining one must 
be chosen to satisfy Equation 8 at s = . 

Next, the general equations of motion in S form will be derived. From Equation 1 5, the 
Lagrangian is 

F = -P - (r 2 + a 2 ) 0 2 sin 2 9 + t 2 - ^ (a 0 sin 2 B + f) 2 . (83) 


Then, from Equation 31, the following four equations corresponding to B, 0, t, and r are obtained: 


d_dF_ _ dF 
ds 30 30 

d_d_ F 3 F 

ds 30 30 

d_ 3F 3 F 

ds 3 t 3 t 

3F _ 3F 
ds dr dr 


= 0, 
= 0 , 
= 0 , 
= 0 . 


(84) 

(85) 

( 86 ) 

(87) 


Equations 84 to 87 and the line-element equation (Equation 15) constitute a set of equations of 
which only four are independent. Equations 15 and 85 to 87 were chosen for this study. In Equation 
83, F is not a function of 0 or t. Hence, Equations 85 and 86 imply that dF/dcp and dF/di are equal to 
constants, denoted respectively as - 2 h and 2e. Therefore, 


— - = - 2 {r 2 + a 2 ) 0 sin 2 0 - 
30 


2 $ra sin 2 0 
P 


( a 0 sin 2 0 + t) = -2 h , 


( 88 ) 



2)3/* 

~F 


{a 0 sin 2 0 + /) = 2e . 


(89) 
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( 90 ) 


These equations are linear in 0 and i and can be solved directly to obtain 

1 h 

0 = [-fiaer + (P - j 3r)] 

PA sin 2 9 


i= —— {(3ahr + e[P(r 2 + a 2 ) + pa 2 r sin 2 6]| , 


From Equation 87, the following Euler equation is obtained corresponding to the r coordinate: 

2P 2 r . r 2 fiQ 2 2 / Br 2 

—r+ — (A? - P A) - — [2rA-P(2r- j3)] - 2r d 2 - 2r 0 2 sin 2 9 - + Q 2 = 0 . (92) 


A A 2 A 2 

Finally, Equation 1 5 differentiated with respect to s becomes 

IP r 2 

— rr+ — (A P- P A) + P d 2 + 2P 9 6 


P P 2 


+ 2r r <p 2 sin 2 9 + 2 (r 2 + a 2 ) 6 <p 2 sin Q cos 9 + 2{r 2 + a 2 ) 0 0 sin 2 6 

2 P P (/ - e) 2 P 

- 2 it + — (/- e) t+ (i- e) 2 r = 0 . (93) 

j3r j 3r fir 2 

Equations 90 to 93 can be put into a form that is convenient for numerical integration by letting 

(94) 

x 2 =e , (95) 

* 3 = 0 , (96) 

* 4 =t, (97) 


Xc=r, 


By definition, we obtain 


x 6 =e 


*1 =*<; , 


( 100 ) 


Then, Equations 90 and 91 become 


*2 = x 6 • 


( 101 ) 


-fiaex 1 + 


sin 2 x. 


(P-pxJ 


( 102 ) 
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1 

x 4 - — — | (3ahx 1 + e[P{x 2 + a 2 ) + Pa 2 x x sin 2 x 2 ] } . (103) 

Solving Equations 92 and 93 for r and 9 and using Equations 94 to 99, it follows that 



x 5 = — 

(104) 


a 21 



a 2i^i~ a n b 2 

X 6 ~ 

(105) 

where 

a l2 a 21 


P = x 2 + a 2 cos 2 x 2 , 


(106) 

A = x 2 - px 1 + a 2 , 


(107) 


5 . 

b 1 = (A P - P A) - Px 2 - 2x l x s <p l sin 2 x 2 - 2(x 2 +a 2 )x 6 <t> 2 sinx 7 cosx 2 

2 P P P 

- 2(x 2 + a 2 ) <j> <j> sin 2 x 2 + 2 it (/- e) t - — (/- e) 2 + (/- e) 2 x 5 , (108) 

jSXj j3x, fix 2 


2x s . _ x 2 p 2 fix 2 

b 2 = (PA - P A) + — [2Xj A - P( 2xj - {3)] + 2XjXg + 2x 1 <j> 2 sin 2 x 2 + — Q 2 - Q 2 , 

A 2 A 2 P P 2 



(109) 

2Px 5 


11 ~ A 

(110) 

= 2Px 
12 ^-*6 ’ 

(111) 

to 

II 

>1*3 

(112) 

P= 2x 1 x 5 - 2a 2 x 6 cosx 2 sinx 2 , 

(113) 

A = 2x 5 x 5 - j3x 5 , 

(114) 

to 

it 

"5? 1 ha 

i 

(115) 
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— { fiahx l + e [P(x 2 + a 2 ) + fia 2 x 1 sin 2 x 2 ] } , 


efiax l + 


sin 2 x. 




( 117 ) 


1 

— | fiahx 5 + e[P(x 2 + a 2 ) + 2x 1 x s P + fia 2 x s sin 2 x 2 

+ 2fia 2 x l x 6 sinx 2 cosx 2 ]| - -^(AP + PA), 


-efiax, + 


(P-frO 


2hx 6 cosx 2 


sin x 0 


(P- fix. ) (PA + PA) . 


(118) 


(119) 


Equivalence of the Constants of Motion in the 5 and E Forms 

In Equation 1 7 for equatorial motion in the E form, there appear the two constants e* and h*, 
which are given by Equations 23 and 24. In the corresponding equation in the S form (Equation 44), 
there appear the constants e and h, which are given by Equations 37 and 38. It is obvious that these 
two sets of constants are not independent of each other. It will be shown that they are indeed identi- 
cal. Moreover, this identity holds not only for the case of motion in the equatorial plane but for all 
motion. 

From the E form of the Kerr metric given by Equation 8, the Lagrangian F*(r *, 0*, r*, 9*, /*, 0*) 


F* = -r* 2 - 2a r* 0* sin 2 9* - ( r * 2 + a 2 ) 0* 2 sin 2 9* 


fir* _ 


- P* 0* 2 + i* 2 - (r* + a 0* sin 2 9* + /*) 2 . (120) 

Because F* is not a function of 0*, it follows that 9P*/90* may be set equal to a constant -2h*\ 


h* = a sin 2 9* 


(l + £) r + sin 2 «* (n« 2 + » 2 + *•+*£-* sin 2 9* . (121) 


Similarly, because F* is not a function of t*, dF*/dt* may be set equal to a constant 2e* to obtain 


fir* . j3 ar* sin 2 9* ■ 




( 122 ) 


Solution of Equations 88 and 89 for motion in S coordinates for h and e yields 

/ n 0 fira 2 sin 2 d\ ■ , n fira sin 2 9 . 

h= I r 2 + a 2 + 1 0 sin 2 9 + t , 
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fi ra ■ ~ 

e = - -jj- 0 sin 2 0 + 


(■-?) 


Next, from Equations 10 to 13, it follows that 


r = r* , 


0 = 0 *, 

0 = 0 * + 

/ = /* - 


A ’ 


/Jr* r* 

A* 


(124) 

(125) 

(126) 

(127) 

(128) 


Substitution for r, 0, 0, and / from Equations 10, 11, 127, and 128 into Equation 123 leads to 

(129) 


, , , |3r*a 2 sin 2 0* 

/z = r* 2 + a 2 + r 




, /. ar*\ (3r*a sin 2 0* / . |3 r* r*\ 

■; sta2 »• r + S + Si — r - — ) 


Rearrangement of terms and substitution from Equation 14 yield 


h = 
or 


a r* sin 2 0 * 
A* 


A* + - + sin 2 0*0* ^r* 2 + a 2 


1 3r*a 2 sin 2 0 *\ fir* a sin 2 0 * . 


+ ; — y, — -) + r ' ' " ' ■ »'* . 030 ) 


, _ •* . 7 , 0r*\ , . 7 j* ( , 7 , sin2 0 *\ . @ r * a sin 2 0* 

h - ar* sin 0 * 11+ — I + sin z 0*0* l r* z + a z 4 I H 

j°* J \ pie j p* 


which is precisely Equation 121 ; hence. 


h — h* . 


(131) 


Similarly, substitution for r, 0, 0, and / from Equations 10, 11, 127, and 128 into Equation 124 
leads to 

fir* r*" 


fir*a 
e = — — - sm 


■ m i e * ( 0 * + ^) + ( 1 - (/* 


p* y a*/ V p*j V A* 

Again, rearrangement of terms and substitution from Equations 14 yield 

j3 r* r* fir* a sin 2 0* . 


e = 


P* 




or 


e- - 


fir* fiar* sin 2 0* . 


which is precisely Equation 122; therefore. 


0 * + 




(132) 


(133) 


— /? ^ 


e = e 


(134) 
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ANALYSIS OF EQUATORIAL ORBITS 


Approach 

A geometric-topological approach is used in the following analysis because there is no known 
direct analytical approach that provides exact solutions. Moreover, even if one exists, it is expected 
that it would be rather long and cumbersome. For instance, in the case of the Schwarzschild gravita- 
tional field, a complete analytical solution of orbits as given by Hagihara (Reference 9) required 1 10 
pages in final journal form and was heavily dependent on elliptic functions. 

As in Appendix A, the results for each value of a are represented in a master diagram that gives 
the complete topological genera of the orbits for various ranges of the two parameters h and e. In 
particular, in the Newtonian limit, the quantity W = e 2 - 1 is the analog of 2 E, where E is the energy 
per unit mass of the test particle. 


Summary of the Equation To Be Studied 

The equation of motion in the S form of a test particle in the equatorial plane of the Kerr gravita- 
tional field has been derived: 


'duS 2 D 2 (u) 

— =B(u) , (44) 

4<t>/ E 2 (u) 


where 

B(u) = j 3(h + ea) 2 u 3 - [h 2 + (1 - e z )a 2 ] u 2 + &u - (1 - e 2 ) , 

D(u ) = a 2 u 2 - (3u + 1 , 

E(u ) = ~P(h + ea)u + h . 

The constants of motion h and e are given by 


h = 



d(j) + fa dt 
ds r ds 


(45) 

(46) 

(47) 


(135) 


e = ■ 


1 3a d(j) 
r ds 


1 


dt 

ds 


(136) 


Without loss of generality, a > 0 is taken for the sake of simplicity of discussion. The quantities 
u, (j), h, and e have in Newtonian orbits the counterparts inverse radius, true anomaly, angular momen- 
tum per unit mass of test particle, and energy per unit mass of test particle, respectively. 


Investigation of the Cubic B(u) 

The cubic B{u) appearing in Equation 44 degenerates into h 2 f(u) for a = 0, where /(w) in this case 
is the one for the Schwarzschild gravitational field and h is a constant. Chan has made a complete 
study of /(;«) for a = 0 and his results are now extended for values of a A 0. In view of this, it is 
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desirable to study the cubic curve B{u). For h + ea^O, B(u) = 0 can have only zero, one, two, or 
three positive roots. A general approach in this direction would require the algebraic knowledge of 
roots of the cubic equation, and this is available, but the analysis becomes very cumbersome algebrai- 
cally. In addition, it becomes more involved, for instance, in the case of one real positive root, 
whether the remaining two roots are real but negative or are complex conjugates. Alternatively, the 
following approach can be taken: On successively differentiating 

B(u) = |3(/z + ea) 2 u 3 - [h 2 + ( 1 - e 2 )a 2 ]u 2 + Bu + e 2 - 1 , (45) 

it follows that 

dB 

B'(u) = — = 3j3(h + ea) 2 u 2 - 2[h 2 + (1 - e 2 )a 2 ] u + j3 ( 137 ) 

du 

d 2 B 

B"{u) = =6P(h + ea) 2 u- 2[h 2 +(1 - e 2 )a 2 ] . (138) 

du 2 

Hence, B"{u) is a straight line intersecting the u axis at 

h 2 +(1 - e 2 )a 2 

u b=u s = , ( 139 ) 

3j 3(h + ea) 2 


where u s is the point of symmetry of the cubic B(u). 

The intercepts, u' Bl and u' B 2 if they exist, of the parabola B'(u) with the u axis are given by 

h 2 + (1 - e 2 )a 2 

t 

U B1,B2 ~ . 

3p(h + ea) 2 

where 

3j3 2 ( h + ea) 2 

V = 1 . (141) 

[h 2 +(1 - e 2 )a 2 ] 2 

Now, in view of equation 139, Equation 140 can also be written in the form 

“'bi,B2=“ s (^VV) 042) 

where for u s > 0, the minus sign is associated with u' B] and the plus sign with u' B2 and for u < 0 the 
plus sign is associated with u R , and the minus sign with u' B2 . 

From the coefficients of the parabola B\u), noting that j3 > 0 in Equation 137, it follows that 
B'(u) can be classified into seven categories, as shown in Figure 1 : 

(a) No real root, u s > 0 

(b) One real positive double root 

(c) Two distinct real positive roots 

(d) No real root, u s = 0 


(i * \/r) 


(140) 
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B'(u) 



Figure 1— Criterion for the seven classes of trace cubic curves. 


(e) No real root, u s < 0 

(f) One real negative double root 

(g) Two distinct real negative roots 

Equations 141 and 142 indicate that 

(1) V< 0 for cases (a), (d), and (e). 

(2) V = 0 for cases (b) and (f). 

(3) 0 < V < 1 for cases (c) and (g). 

Corresponding to these, the cubic B(u) can be classified into seven categories as shown in Figure 
2; for each case, all conceivable curves are qualitatively illustrated and are not obtained by simple 
vertical translation from one another as can be verified by Equation 45. As a concrete example for 
discussion, the cubic-curve category c was chosen because it contains many features shared by the 
other types. The cubic B(u) is shown in Figure 3 with the three roots u BX B2 53 (indicated by w, n m 
in the figure) of B(u) = 0. In general, any real root of B(u) = 0 is denoted by u B . 

For a double root, 

B(u b )= 0 (143) 


B'(u b ) = 0. 

Hence, double roots of two types occur: 


( 1 ) 


U B 1 U B2 U B1 ' 


(144) 

(145) 
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B(u) B(u) 




Figure 2— Curves for each of the seven classes of trace cubic. 
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Figure 3— Cubic curve of category c with W< 0. 


(2) u B2 Wg 3 u B2 . (146) 

For double roots of type 1 ,u B1 = u B2 and Equation 143 becomes B(u' bi ) = 0. Substituting u' Bl 

from Equation 142, 


- s/V J^- (3u s - 2(3(h + ea) 2 u\ 

the square of which is 

"2 


( 1 - c 2 ) - j 3u + 2j3 {h + ea ) 2 u ~ , 


V 


3 j3u s - 2 j5(h +ea) 2 u s 3 


= [(1 ~ e 2 ) - fiu + 2(3(h + ea) 2 u 2 ] 2 . 


(147) 


(148) 


Similarly, for double roots of type 2, u B2 = u' B2 , and Equation 143 becomes B(u' B2 ) = 0. Substi- 
tuting u' S2 from Equation 142 and then squaring the resultant equation yields again Equation 148. 
That is, Equation 148 gives solutions of both type 1 and type 2 double roots of B(u) = 0. 

Substituting for u s and V from Equations 139 and 141 into Equation 148, the following equa- 
tion is obtained after much tedious algebraic work: 

A 6 h 6 +A s h 5 +A 4 h 4 +A 3 h 3 +A 2 h 2 +A x h+A Q = 0 , (149) 

where the coefficients A f (e, a, j3) are given by 

A Q = a 2 [4(1 - e 2 )V +/3 2 a 2 (27e 8 - 72c 6 + 62c 4 - 16c 2 - l) + 4/3 4 e 2 ] , (150a) 

A 1 =P 2 ea[a 2 (l08e 6 - 252e 4 + 180c 2 - 36) + 8/3 2 ] , 


(150b) 



(150c) 


A 2 = 12(1 - e 2 )V + p 2 a 2 ( 162e 6 - 324e 4 + 182e 2 - 20) + 4 p 4 , 

A 3 = j3 2 ea(108e 4 - 1 80e 2 + 72) , (150d) 

A a = 12(1 - e 2 ) 2 a 2 + (3 2 (8 - 36e 2 +27e 4 ) , (150e) 


^5 = 0 ’ 


(1500 


A 6 =4(1 - e 2 ). 


(15 Og) 


Now, in principle, h can be solved as a function of e with a and |6 as parameters. However, this 
cannot be done in general because Equation 149 is of the sixth order. Instead, a relation between h 
and e with a and j3 as parameters may be found by numerical methods. That is, for a given set ( a , /3) a 
value of e is fixed and then values of h that satisfy Equation 149 are computed. This is then done for 
another value of e to obtain the corresponding values of h. Some typical results are presented in 
Figures 4 to 7; these figures are for |3 = 2 and a = 0, 0.7, 1 , and 1 .4, respectively. The symbols S and U 
designate, respectively, the loci for which double roots of 5(«) = 0 of types 1 and 2 occur. That is, the 
S line is given by eliminating u' B ] between 


and 


B(u' b1 ; h, e, a, J3) = 0 
B'(u' Bl ;h,e,a,P) = 0 , 


(151) 


whereas the U line is given by eliminating u’ B2 between 


B(u' B2 ;h,e,a,P) = 0 

and 


(152) 


B'(u' B2 ;h,e,a,@) = 0 . | 

It follows from Equations 142, 145, and 146 that there are no double roots of B(u ) = 0 for 
V < 0. Because of this, the curves for which V = 0, known as V lines, are also plotted in Figures 4 to 
7. From Equation 45 it can be seen that the intercept of the cubic B(u) with the B(u) axis is given by 


5(0) = <? 2 - 1 = W , 

so that 

5(0) <0 for W < 0 

5(0) >0 for > 0 . 


In these figures, the W lines are defined by W(e) = 0, where W(e) is given by Equation 153. 


(153) 

(154) 


Next, the cubic B(u) is classified according to regions of the h, e space for a given set of (a, j3). 
For a = 0, B(u) degenerates to h 2 f(u ) where /(w) is the cubic studied by Chan in Appendix A. In 
Figure A15 and Table A2, he has already classified the cubic into seven essentially different regions of 
the h, W space labeled 1, 2, 3, 4, 4*, 5, and 5*. Table 1 gives the properties of the cubic B(u) In these 
regions and two other regions 6 and 6* which exist because of the additional cases of -the quadratic 
B’(u) corresponding to u s < 0. 
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Degeneracy of V, U, U, V Lines 
for All Values of e and Degeneracy 
of Z Lines for e 2 Sl. 



SSdW ilNfl y3d A9SJ3N3 
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Figure 4-Classification of h, e space according to characteristics of B(u) curve for a = 0.0. 
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Table 1 —Characteristics of the cubic B(u ) 
for various regions of the master diagram. 


Region 

B ( u \) 


V 

w 

u s 

1 

<0 

<0 

>0 

<0 

X 

2 

>0 

<0 

>0 

<0 

X 

3 

>0 

<0 

>0 

X ) 

X 

4 

>0 

X ) 

>0 

>0 

X 

4 * 

C) 

( a ) 

<0 

>0 

X 

5 

>0 

>0 

>0 

<0 

X 

5 * 

( a ) 

( a ) 

<0 

<0 

X 

6 

>0 

>0 

>0 

X ) 

<0 

6* 

( a ) 

( a ) 

<0 

X ) 

<0 


and u'g 2 do not exist for V < 0. 


These cubics are shown in Figure 8. Using a one-to-one correspondence between Figures 4 and 
A15, the regions in Figure 4 are labeled. A little consideration reveals that the regions in Figures 5 to 
7 can be labeled as shown. Figure 4 is a master diagram for a = 0; however, Figures 5 to 7 are not yet 
master diagrams because f(u) = B(u)D 2 (u)/E 2 (u) and not just B(u)/h 2 must be considered fora 0. 

For convenience, Z and H are defined as 

Z = h 2 +(1 - e 2 )a 2 (155) 

H = h + ea . (156) 

For Z > 0 and H^O, the cubic B(u) of Equation 45 is similar to the one studied by Chan, which is 

2m W 

f(u ) = 2 mu 3 - u 2 + — u + — (157) 

h 2 h 2 

and hence the labeling of regions 1, 2, 3, 4, 4*, 5, and 5* of Figures 5 to 7 follow from Chan’s work 
(Appendix A). The characteristics of the cubic B(u) corresponding to these regions are shown in 
Figure 8, which is equivalent to Table 1. We have already noted that for all seven of these regions the 
point of symmetry u s given by Equation 139 and the roots u' Bl and u' B2 of B'(u) = 0 given by Equa- 
tion 140 are greater than zero. 

In addition, two new regions previously mentioned occur for Z < 0 and H =£ 0. In these two 
regions the point of symmetry of the B(u ) curve and the roots of B'(u) = 0 are both less than zero. In 
region 6*, since V < 0, it is obvious that B(u) = 0 has only a single root, which is negative. In region 6, 
it can be proved that B ( u ' b1 ) > 0 and B ( u B2 ) > 0 so that there can only be one negative root of 
B(u ) = 0. The details are given in Appendix B. Because of these additional regions, the lines for which 

Z = h 2 + ( 1 - e 2 )a 2 = 0 , (158) 

known as the Z lines, are plotted in Figures 5 to 7. 
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Figure 8— Characteristics of the B(u) curve in various regions. 


Finally, it remains to discuss the special case in which H = 0. For this, it is convenient to define 
the H line by 

H = h + ea = 0, (159) 

which is always a straight line. From Equations 45 and 159 it is seen that on an# line, B(u) becomes 
the quadratic 

B{u) = -a 2 u 2 + $u + e 2 - 1 . (160) 

The roots of this equation are 

(161) 
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From Equation 161, it is noted that Equation 160 has distinct roots for 

P 2 

e 2 >\ , (162) 

4a 2 

a double root at 

P 2 

e 2 = 1 , (163) 

4a 2 

and no real roots for 

P 2 

e 2 <1 . (164) 

4a 2 

As a summary, in the h, e plane for a given set ( a , p)— 

( 1 ) The S and U lines give the condition for double roots of the cubic B(u) = 0. 

(2) The V lines give the condition for double roots of the quadratic B'(u) = 0. 

(3) The W lines give the condition for B( 0) = 0. 

(4) The Z lines give the condition for which the point of symmetry u g = 0. 

(5) The H lines give the condition for which B(u ) is a quadratic; otherwise, B{u) is a cubic. 


Investigation of the Function f(u) 

Now consider the equation 

D 2 (u) 

m=B(u ) ( 165 ) 

E 2 {u) 

where B(u), D(u), and E(u) (defined in Equations 45 to 47 ) are cubic, quadratic, and linear in u, 
respectively. Let Uj be any root of the equation 

/( h ) = 0 . ( 166 ) 

For a double root, 

f(u f ) = 0 ( 167 ) 

f\u f ) = 0 . ( 168 ) 

Specifically, let u £ be any value for which these two equations hold. It can be shown that Equations 
167 and 168 are the necessary and sufficient condition for a circular orbit to exist. (See Appendix A 
discussion of Equations A7 to A10 and A 19.) 

Table 2 presents all of the conceivable combinations of the zeros of B(u), D(u), and E(u) for 
which Equation 166 could possibly have a degenerate root u £ . That most of these cases are impossible 
can be shown as follows: Let u B ,u D , and u £ be any real root of B(u) = 0, D(u) = 0, and E{u) = 0, 
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Table 2— Ways in which f(u ) = 0 could have a degenerate root u c . 


Case 

Number of Roots of 
B(u) = 0 Equal to u c 

Number of Roots of 
D(u ) = 0 Equal to u c 

Number of Roots of 
£■(«) = 0 Equal to u c 

l a 

0 

1 

0 

2 a 

0 

2 

0 

3 

0 

2 

1 

4 

1 

1 

0 

5 

1 

2 

0 

6 

1 

2 

1 

7 a 

2 

0 

0 

8 

2 

1 

0 

9 

2 

1 

1 

10 

2 

2 

0 

ll a 

2 

2 

1 

12 a 

3 

0 

0 

13 

3 

1 

0 

14 

3 

1 

1 

15 

3 

2 

0 

16 a 

3 

2 

1 


a These cases axe actually possible. 


respectively. That is, by definition 


The condition 


B(u b ) = 0 , 
D(u d ) = Q , 
E(u e ) = 0. 
B(u d ) = 0 


U B U D 


is examined first. From Equation 49, which is 


AD + a 2 B=E 2 , 

and the definition D(u d ) = 0, 

a 2 B(u D ) = E 2 (u D ). 

It follows that 

B{u D ) = 0 <=> E(u d ) = 0 . 
Hence, cases 3, 4, 5, 8, 10, 13, and 15 of Table 2 are impossible. 
Next, the condition 


( 169 ) 


( 170 ) 


( 49 ) 


(171) 

( 172 ) 


( 173 ) 


D'(^) = 0 
B(u D ) = 0 
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is considered. Differentiation of Equation 49 results in 

A'D+AD' + a 2 B' = 2EE ' . 

From Equations 49, 174, and the definition D(u d ) = 0, it follows that 

A(u d ) D\u d ) + a 2 B'{u d ) = 2E(u d ) E'(u D ) . 
Hence, from Equations 171, 173, and 175, 

D'(u d ) = 0 


=*B'{u d ) = 0. 


B(u d ) = 0 

Therefore, case 6 of Table 2 is impossible. 

The condition 


( 174 ) 


(175) 


(176) 


B(u £ ) = 0 
B'(u e ) = 0 


(177) 


is examined next. From Equations 49 and 174 and the definition E(u e ) = 0, the following are obtained: 

A(u E )D(u E ) + a 2 B(u E ) = 0 (178) 

A'(u e )D(u e ) +A(u e ) D'(u e ) + a 2 B'(u e ) = 0 . (179) 

From Equations 177 to 179, therefore, 

A(u e )D(u e ) = 0 (180) 

A\u e )D(u e ) +A(u e ) D\u e ) = 0 . (181) 

For Equation 180, there are three cases: (a) A{u E ) = 0 and D(u E ) = 0, (b) A(u E ) = 0 and 
D(u e ) =£ 0, and (c) A(u E ) 0 and D(u E ) = 0. It is shown in Appendix B that A (u E ) = 0 and 
D(u £ ) = 0 cannot be satisfied simultaneously and hence case a is impossible. From Equation 181, it 

is seen that 

A\u e )D{u e ) = 0, (182) 

so that 

A'(u e )= 0. (183) 

It is shown in Appendix B that A'(u £ ) ¥= 0 for any finite u E , and hence case b is impossible. Finally, 

from Equation 181, 

< 184 ) 

so that 


It follows from this discussion that 


A(u e )D'(u e ) = 0 
D\u e ) = 0. 


B(u e ) = 0 
B'{u e ) = 0 


D(u e ) = 0 
( D'{u e ) = 0 


(185) 


( 186 ) 



From the identity 
Equation 186 becomes 



(187) 

(188) 


Hence, cases 9 and 14 of Table 2 are impossible. 


From Equations 45 to 47, we see that B(u;h, e, a, j3) = 0, D(u;a, (3) = 0, and E(u; h, e, a, /?) = 0 do 
not contain any common factors for all values of ( h , e, a , (3). It follows that there exist values of 
( h , e, a, (3) for which cases 1, 2, 7, and 12 are possible. 


From Equations 172 and 176, the following is obtained: 

D(u d ) = 0 

D'(u d ) = 0 


=> < 


B(u d ) = 0 
B'{u d ) = 0 


(189) 


E(u d ) = 0 


Because it is obvious from Equations 46 and 47 that the left-hand side of Equation 1 89 can be satisfied, 
case 1 1 is possible. 


Finally, for case 16, consider only the equations 


D{u d ) = 0 , 
D\u d ) = 0, 
E(u d ) = 0 , 
B"(u d ) = 0. 


(190) 


The remaining two equations, B(u D ) = 0 and B'(u D ) = 0, are automatically satisfied in view of Equa- 
tion 1 89. It is shown in Appendix B that this set of equations has the solution given by the following 
equations: 


1 


e = ± x/3 ’ 

(191a) 

2 

+ % /3 a ’ 

(191b) 

= ±2a , 

(191c) 

1 

* D ± a ’ 

( 1 9 1 d) 
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Hence, case 16 is possible. 

In summary, all the possible cases for which /(u) = 0 can have a degenerate root at u c are 1, 2, 7, 
11, 12, and 16, as indicated in Table 2. 


Next, the properties of the cubic B(u ) at the values u D and u E will be investigated. From Equa- 
tion 171 it follows that 


1 


B(u b )‘-EHu d ) 

a 1 


[h - PQi + ea)u D ] 2 


1 

= — Qi ~ h$u D - a@eu D ) 2 , 
a 2 

which, in view of the definition D(u d ) = 0 or 

D(u d ) = a 2 uj ) - (Su D + 1 = 0, 


becomes 


1 


It is obvious that 


B(u d ) = — [h - h{a l Up + 1) - afieu D ] ' 
a 2 

= u 2 d ( ahu D + &e) 2 . 

e/3 

B(u n ) = 0<=*h = - — 

nil 


au 


D 


B(u d ) > 0 otherwise , 


(192) 


(193) 


(194) 


(195) 


where u D is either of the two real roots u D l and u D 2 . 

In view of this, the D { line and D 2 line in the h, e space for a given set (a, j3) are defined by 

the respective equations: . 

ep 


D { = h + 


D 2 = h + 


au 


D 1 


eP 


au 


= 0 


= 0 , 


D 2 


(196) 


where u D1 and u no are the two roots of D(u) = 0; i.e., 


D2 


P ± sjp 2 - 4 a 2 


u 


D1,D2 


(197) 


2 a 2 
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In general, a D line refers to either one of these and thus may be defined by eliminating u between the 
equations 




B(u;h, e, a, P) = 0 
D(u;a, p) = 0 . 

It is shown in Appendix B that B(u F ) has the following two equivalent forms: 


ae~ 


B(u E ) = 


B(u f ) = 


P 2 (h + e) 2 


C eh - a + ae 2 ) 


(h\ 2 j3 2 h 

+ p 2 

\ej a e 


(eh - a + ae 2 ) (h + 


ej3 


h + 


ep 


au 


D 1 / 


au 


Dl 


(198) 


(199) 


( 200 ) 


P 2 ( h + ae) 2 

For convenience, the E x line in the h, e space for a given value of a is defined as 

E x = eh - a + ae 2 = 0 . 

For the case p 2 - 4a 2 > 0, so that by Equation 197 the roots of D(u) = 0 are real, it is preferable 
to use Equation 200 from which results 


( 201 ) 


B(u e ) 0 


h, e on D line 
^ or h,e on E x line 


( 202 ) 


B(u e ) ^ 0 otherwise. 


In this case, the sign of B(u E ) is determined by the combination of the D and E x lines. Typical results 
are presented in Figures 9 and 10 for P = 2 and a = 0.7 and a = 1, respectively. In Figure 10 the regions 
between the D lines disappear for a = 1. 

For the case p 2 - 4 a 2 < 0, so that the roots of D(u) = 0 are not real, it is preferable to use Equa- 
tion 199 from which can be seen that 


in view of 


because the discriminant (p 2 /a 2 )(p 2 


Hence, in this case the sign of B(u E ) 
in Figure 11 for P = 2 and a = 1.4. 


B(u e ) = 0 ■*=> h, e on E x line 


(203) 


0 ) 


2 0 2 h , 

+ + p 2 0 


a e 

4a 2 ) <0. Moreover, 


Q 


h\ 2 P 2 h , 

-1 + — - + p 2 > 0 . 
a e 


(204) 


is determined by the E x line only. Typical results are presented 


31 








B(u e )>0 
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and 

but 


It follows that the £’ 1 line is defined by eliminating u between 

B(u,h, e, a,(S) = 0 

E(u;h,e,a,i 3) = 0 ; 

D(u ;a, (3) =£0 . 


(205) 


Now that the sign of B(u D ) and B(u E ) can be determined in various regions of the h, e space, the 
conditions that determine the relation between the order of u £ ,u D1 , and u D 2 must be found as follows: 


From the definition of u E which yields 

it is noted that the difference 5 defined by 
can also be written as 


h 


U E 


(3(h + ea ) 


8 —u D u E 


8 = 


1 


/3(/z + ea) 


[$u D ea + h{$u D - 1)] 


where u D is defined by 

a 2 u E - ($u D + 1=0 

From Equations 208 and 209, it follows that 

/Se' 


5 = 


a 2u l 


\h + 


(3(h + ea) 


au T 


or 


2„2 


an 


~ U D 1 U E 


2)1 


)3(/z + ea) 


2„2 


a“u 


^2 ~ U D2 U E 


D 2 


\h + 


h + 


$(Ji + ea) 



Thus, the ordering of u E ,u D1 , and u D2 is determined by the H and D lines. 
Moreover, from Equation 209 and the relation 


U Dl U D2 


1 

n2 


(206) 

(207) 

(208) 

(209) 

( 210 ) 


( 211 ) 


( 212 ) 


(213) 


derived from it, Equation 196 for the two D lines may also be written as 

Z> 1 = h + ea(l + a 2 u E2 ) = 0 | 

D 2 =h + ea(l + a 2 u El ) = 0 . ) 

Hence, the H line does not lie between the two D lines. Also, the Z> 1 line has a less-negative slope than 
th qD 2 line. Typical results are presented in Figures 12to 14for (S= 2 anda = 0.7, l,and 1.4, respectively. 


35 









Figure 14-Typical classification of h, e space according to the ordering of u £ , u D1 , and u D2 and sign of u E for 1 < a < °° (shown here fora = 1.4); 

in the case 1 .0 < a < °°, no roots o D1 and u D2 can exist for D(u) = 0. 



(In these diagrams, subscripts 1 and 2 in the two D lines and subscript 1 from the E l line have 
been omitted for the sake of simplicity. This use of notation should not result in the erroneous 
interpretation of these lines as the D(u ) = 0 and E(u) = 0 lines, as defined in Equations 46 and 47. 
Rather, any line in the h, e space can denote at most only a relation F(h, e; a, (3) = 0, where a and jS 
are parameters, and cannot involve u.) 


Finally, from Equations 197 and 206 it can be seen that 



U D1,D2 

>0, 


>0 

for 

h 

>0, 

h + ea 

<0 

for 

h 

<0, 

h + ea 

= 0 

for 

h 

= 0, 

= +00 

for 

h + ea 

= 0 . 


(214) 


(215) 


As a summary of this investigation of the function /(«), the following points are given concerning 
the h, e space for a given set ( a , (3): 

(1) The positive definiteness of B(u D ) is determined by the D lines. 

(2) The sign of B(u £ ) is determined by the combination of the D and E l lines. 

(3) The ordering of u £ ,u D1 , and u D2 is determined by the combination of the D and H lines. 

(4) u Dl and u D2 depend only on (a, {3); u D1 >0 and u D2 >0. 

(5) The sign of u E is determined by the combination of the e axis, h = 0, and the H line, 
h + ea = 0. 

(6) There are two D lines for (3 2 - 4a 2 > 0, one degenerate D line for j3 2 - 4a 2 = 0, and none for 
P 2 - 4a 2 <0. 


Construction of the Master Diagram 

In this section, the results of the two preceding sections will be used to construct the master dia- 
gram for each a such that it will give all admissible orbits corresponding to a given set h, e and /3 - 2 in 
Equation 44. 

For concreteness, first consider the case of a = 0.7, which illustrates all the features of 0 < a < 1. 
Superposition of Figures 5, 9, and 12 produces Figure 15. It is now necessary to divide the regions 1, 
2, 3, 4, 4*, 5, 5*, 6, and 6* into subregions which we will label with a letter as a suffix (Table 3). With 
the information in Figure 8 and Table 3, it is quite easy to determine the ordering of the roots of 
B(u) = 0, D(u) = 0, and E(u) = 0 in the subregions of 1, 4, 4*, 5, 5*, 6, and 6*. However, this deter- 
mination is more difficult in the subregions of 2 and 3 because it appears that the values of all these 
roots must be known explicitly to effect their ordering. The roots of D(u) = 0 and E(u) = 0 are given 
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Table 3 -Characteristics of the cubic B(u) for various subregions of h,e space for 0 < a < 1 . 


Region 

U E 

Ordering of 
U D1’ U D2’ U E 

B(ug) 


B(u' n ) 

V 

w 

U S 

la 


>0 

C4 

a Q 

V 

a Q 

V 

£ 

>0 





X) 

lb 


>0 

U E <U D1 <U D2 

<0 

<o 

<0 

>0 

<0 

lc 


>0 

u Dl <u E < - u D2 

>0 






2a 


>0 

U E < ' u Dl < ' u D2 

>0 






2b 


>0 

U D\ < U D2 < ' U E 

<0 






2b' 


<0 

U E <u Dl < ~ u D2 

<0 

>0 

<0 

>0 

<0 

X) 

2c 


>0 

U D\ <U E < ' U D2 

>0 




2d 


>0 

U D1 < ' U E < ' U D2 

>0 






2e 


>0 

<N 

£ 

V 

T“*i 

£ 

V 

a * 

<0 






3a 


>0 

u E <u Dl <u D1 

>0 






3b 


>0 

U D 1 <U D2 < ' U E 

<0 






3b' 


<0 

U E < ' U D1 < ' U D2 

<0 

>0 

<0 

>0 

X) 

X) 

3c 


>0 

U m < u £ < U D2 

>0 



3d 


>0 

u Dl < ' u E <U D2 

>0 






3e 


>0 

U E < ' U D1 < ' U D2 

<0 






4a 


>0 

cq 

I 5 

V 
*—■ 

V 
£ 

>0 






4b 


<0 

U E < ~ U D2 

<0 

>0 

>0 

>0 

X) 

X) 

4c 


>0 

u Dl < U E < ' U D2 

>0 






4a* 


>0 

U E < ' U D1 < ' U D2 

>0 






4b* 


<0 

U E < ' U D\ < ' U D2 

<0 

( a ) 

( a ) 

<0 

X) 

X3 

4c* 


>0 

u Dl <U E < U D2 

>0 






5a 


>0 

U E < ' u Dl < - u D2 

>0 






5b 


<0 

U E < - u Dl < ' u D2 

<0 

>0 

>0 

X) 

<0 

X ) 

5c 


>0 

U E < ' u Dl < ' u D2 

>0 






5a* 


>0 

U E < ' U D1 < ' U D2 

>0 






5b* 


>0 

U E < - u Dl <u D2 

<0 

( a ) 

( a ) 

<0 

<0 

X ) 

5b*' 


<0 

U E < u Dl < ' u D2 

<0 

5c* 


>0 

U D\ < - U E <U D2 

>0 






6a 


>0 

cq 

S 

V 

r— ( 

:? 

V 

a * 

>0 

>0 

>0 

>0 

X ) 

<0 

6b 


<0 

U E <U D1 <U D2 


6a* 


>0 

C4 

a Q 

V 

V 
a * 

>0 

( a ) 

( a ) 

<0 

X ) 

<0 

6b* 


<0 

U E < U D1 <U D2 

so 


a Because V < 0 in these regions, u'j arid i/jj do not exist. 

Note: Because in subregions 2b, 2b', 3b, 3b', 5b*, and 5b*' B(u E ) < 0 and (du/d<p) 2 = f\u E ) = B(u e )[D 2 (u e )/E 2 (u e )] , it is 
unnecessary to distinguish between the primed and unprimed subregions in the master diagram. 


explicitly by equations 197 and 206. The roots of the cubic B{u) = 0 are available from the algebraic 
theory of cubic equations but, as already mentioned, this leads to very cumbersome algebraic expres- 
sions. Alternately, the same objective may be achieved by going continuously along some desirable 
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curve in the h, e space from a subregion with known ordering to one whose ordering has to be 

determined. 

For instance, consider the boxed region P in the third quadrant of Figure 1 5. If we begin with 

(a) any point in subregion la, then we know from Figure 8 and Table 3 that the ordering of the roots 
is as shown in Figure 16a. Let us now take any continuous curve that goes from this starting point to 

(b) any point on the boundary given by the S line, then to (c) any point in region 2a, next to (d) any 
point on the boundary given by the W line, and finally to (e) any point in subregion 3a. Then, the 
ordering of the roots is obtained as shown in Figure 16. 





B(u) 

/"\ /, , 

i 

7 

V y ' ' 

u i u n v — u m u e u m 

1 

U D2 


(d) W-LiNE BETWEEN SUBREGIONS 2a AND 3a 



Figure 16— Ordering of u E , u D1 , and u D2 in some 
illustrative subregions of the third quadrant in the 
master diagram of Figure 15. 


In the same way, we obtain the ordering in 
subregions 2d, 2e, 3d, and 3e in the fourth quad- 
rant of Figure 1 5 by starting from any point in 
subregion lc as shown in Figure 17. 

Next, to determine the ordering of the roots 
in subregions 2b and 2c in the second quadrant 
of Figure 15, it is simplest to consider first any 
point on the H line in subregion 2b. The Bin) 
curve is now the degenerate quadratic given by 
Equation 160. From Equation 195, B(u D ) > 0 
because the point under consideration is not on 
a D line. From Equation 2 1 5, it follows that 
u E = ±°°. Consequently, the ordering of the 
roots is as shown in Figure 1 8a. As described for 
Figure 16, the ordering of the roots is obtained 
as shown in Figure 1 8. 

Finally, the ordering of the roots in subregions 
3b and 3c can be obtained from Figure 18; the 
ordering is the same except that now u Bl < 0 be- 
cause J3(0) > 0 in region 3. 

From a knowledge of the ordering of the roots 
of B{u ) = 0, D(u ) = 0, and E{u ) = 0 and from the 
fact that 

D 2 {u) 

m=Biu) , (165) 

E 2 (u) 

a complete characterization can be made of the 
f(u ) curve in various subregions of the h, e space 
for a given set (a, j3). For the present case of 
0 < a < 1, this characterization is given in Fig- 
ure 19. For this reason, Figure 15 is known as a 
“master diagram.” It is specifically for a = 0.7, 
but qualitatively valid for 0 < a < 1. 
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(a) SUBREGION lc (b) S-LINE BETWEEN SUBREGIONS (c) SUBREGION 2d 

1c AND 2d 



(g) SUBREGION 3e (h) W-LINE BETWEEN SUBREGIONS (0 SUBREGION 2e 

3e AND 2e 


Figure 17— Ordering of u E , u 01 , and u D2 in some illustrative subregions of the fourth quadrant in the master diagram 

of Figure 15. 



(a) H-LINE BETWEEN SUBREGIONS (b) SUBREGIONS 2b (left of H line) (c) SUBREGION 2b' (right of H line) 

2b AND 2b' 



(d) D-LINE BETWEEN SUBREGIONS (e) SUBREGION 2c 

2b AND 2c 

Figure 18— Ordering of u E , w 01 , and u D2 in some illustrative subregions of the second quadrant in the master diagram 

of Figure 15. 
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Figure 19 (continued)— Characteristics of the f(u) curve in various subregions for 0<a < 1. 
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REGION 3b- 3c (D2 line) REGION la-lb (E line) 



f(u) CURVE f(u) CURVE 

REGION 2b (H line) REGION 3b (H line) 



REGION 1c-2d IS line) REGION 2a-5a (U line) 


"'-3 


Figure 19 (continued)— Characteristics of 



REGION 6b* (E line) REGION 5a*-5b* (E line) 



REGION la- 2a IS line) REGION 1b-2e IS line) 



f(u) curve in various subregions for 0 < a < 1 , 
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In the same way, the characterization for the case of a = 1 is given in Figures 20 and 21 and 
Table 4, and for the case a= 1.4, which is qualitatively valid for 1 <a < the characterization is 
given in Figures 22 and 23 and Table 5. 

Master diagrams for |3 = 2 and for values of a ranging from a = 0 to a - 100 are presented in Appen- 
dix C. Figures Cl to CIO are for a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9; Figure Cl 1 is for 
a =1; and Figures C12 to C38 are fora = 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 3, 4, 5, 6, 7, 8, 9, 
10, 20, 30, 40, 50, 60, 70, 80, 90, and 100. 


The method of obtaining the mathematical solutions u(<p) from the differential equation, Equa- 
tion 44, is discussed next. The equation 

/du\ 2 D 2 ( u ) 

— =m=B{u ) (216) 

\d(p/ E 2 (u) 

is autonomous; therefore, the whole one-parameter family of solutions is given by 

F(<p + c,u) = 0, (217) 

where c is an arbitrary constant. That is, if any one solution in the family is known, the others in the 
family are obtained by a translation. However, from the theory of nonlinear first-order ordinary dif- 
ferential equations it is possible that there are singular solutions of Equation 216 that are not included 
in Equation 217. These singular solutions are of various types, among which are the “isolated” solu- 
tion, the envelope, and the “asymptotic” envelope as discussed in Appendix A. Hence, because the 
f(u) curve is already classified in Figures 19, 21, and 23 for 0<a< 1 , a = 1, and 1 < a < °°, respec- 
tively, it is quite easy to sketch qualitatively the solutions u(<P) in the various subregions of the three 
kinds of master diagrams typified by Figures 1 5, 20, and 22. 


The special case in which there is a subregion with f(u E ) = +°°; i.e., B(u E ) > 0, D(u E ) 0, and 

E(u e ) = 0, gives quite unexpected results because it leads to 



(218) 


and, at first sight, it is uncertain which of the two branches should be taken to continue the solution 
u(0) at u E . The governing principle is that the solution is to be continued in such a way that the proper 
time s always increases. From Equations 14, 42, 46, and 47, the following is obtained: 


d(j) _ u 2 E{u ) 
ds D{u) 


(219) 


From this equation and the condition given by Equation 218, it follows that an orbit must be chosen 
with a reversal in <j> at u E . 

The qualitative solutions for each subregion of the master diagram are presented in Appendix D 
for 0 <a < 1 (Figures D1 to D15), a = 1 (Figures D 16 to D29), and 1 <a <°° (Figures D30 to D42). 
In addition, the “physical orbits”, i.e., orbits in the r, <j> plane, are sketched. 
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Figure 20-Master diagram for classification of orbits for a = 1.0, obtained from superposition of Figures 6, 10, and 13. 








REGION 2e REGION 3a 




REGION 3b REGION 3e 



REGION 4a* REGION 4b* 


curve in various subregions fora = 1. 
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Figure 21 (continued)— Characteristics of the f(u) curve in various subregions for a = 



B(u) 
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Figure 21 (concluded)— Characteristics of the f{u) curve in various subregions for a = 



Table 4— Characteristics of the cubic B(u) for various subregions of h, e space for a = 1. 


Region 

U E 

Ordering of 

U D =u Dl ~ U D2 andz/j. 

B( u E ) 


5(nJj) 

V 

w 

u s 


>0 

u E <u D 

X) 







>0 

u E <u D 

<0 

<0 

<0 

X3 

<0 

X) 










2a 

>0 

U E <U D 

X) 






2b 

X) 

u D <u E 







— 

<0 

u E <u D 

<0 

X) 

<0 

X) 

<0 

X) 

mm 









2e 

>0 

U E ^ U D 





‘ 


3a 

>0 

U E < ' U D 

X) 






3b 

>0 

U D < ' U E 







3b' 

3c 

<0 

U e <u d 

<0 

X) 

<0 

X) 

X) 

X) 

3d 









3e 

>0 

U E <U D 

<0 






4a 

>0 

U E < ' u D 

X) 






4b 

<0 

U E < - u D 

<0 

X) 

X) 

X) 

X) 

X) 

4c 









4a* 

>0 

u E <u D 

X) 






4b* 

<0 

U E < U D 

<0 

( a ) 

( a ) 

<0 

>0 

X) 

4c* 









5a 

>0 

U E < ' U D 

X) 






5b 

<0 

U E <U D 


X) 

X) 

X) 

<0 

X) 

5c 



<0 






5a* 

>0 

U E < U D 

X) 






5b* 

>0 

U E < ' U D 

<0 

( a ) 

( a ) 

<0 

<0 

X) 

5b*' 

<0 

U D < ' U E 







5c* 









6a 

>0 

U E < - u D 

X) 

X) 

X) 

X) 

X) 

<0 

6b 

<0 

U E < U D 

X) 






6a* 

X) 

U E < ' U D 

X) 

( a ) 

( a ) 

<0 

X) 

<0 

6b* 

<0 

U E < ' U D 

X) 







a Because V < 0 in these regions, uj and Ujj do not exist. 

Note: Because in subregions 2b, 2b', 3b, 3b', 5b*, and 5b*' B(ug) < 0 and ( du/d<p ) 2 = /(«£•) = B(u e )[D 2 (u e )IE 2 (u e )], it is 
unnecessary to distinguish between the primed and unprimed subregions in the master diagram. 
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Figure 22-Typical master diagram for classification of 








f(u) CURVE 


REGION 3a REGION 3b 



f(u) CURVE 
REGION 4a 



REGION 4a* REGION 4b* REGION 5a 




REGION 5b 


Figure 23— Characteristics of the f(u) curve in various subregions for 1 < a < 
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Table 5— Characteristics of the cubic B{u) for various subregions of 

h, e space for a > 1 . a 



a Because a > 1 and p = 2, u Dl and u D2 do not exist. 
b Because V < 0 in these regions, uj and ujj do not exist. 

Note: Because in subregions 2b, 2b', 3b, 3b', 5b*, and 5b*' B(u E ) < 0 and 
( du/d(t > ) 2 =/(«g) = B(u e )[D 2 (u e )IE 2 (u e )] , it is unnecessary to distinguish between the primed 
and unprimed subregions in the master diagram. 
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From these diagrams, a collection of all the various types of orbits is obtained (Figure 24); there 
are 25 essentially different ones. In this diagram, the name given to describe each orbit attempts to 
illustrate its properties at r max and r mjn . For example, for orbit of type 6, the outer part of the orbit 
is like that of a parabola whereas the inner part is that of a limit circle. 

All the results of the solutions given in Appendix D and Figure 24 are tabulated in Table 6, which 
shows all admissible orbits in each subregion of a master diagram for 0 < a < 1 , a = 1 , and 1 < a < °°. 



TYPE 6: PARABOLIC LIMIT CIRCLE 

TYPE 5r: ELLIPTIC LIMIT CIRCLE 
WITH REVERSAL 




TYPE 4: HYPERBOLIC -LIKE ORBIT 



TYPE 5: ELLIPTIC LIMIT CIRCLE 



PAKABULIU LiM! 
WITH REVERSAL 



TYPE 8: LIMIT-CIRCLE SPIRAL 



TYPE 10r : PARABOLIC SPIRAL 


WITH REVERSAL 



TYPE 8r: LIMIT-CIRCLE SPIRAL 
WITH REVERSAL 



TYPE 11: HYPERBOLIC SPIRAL TYPE lir: HYPERBOLIC SPIRAL 


WITH REVERSAL 



TYPE 13: LIMIT-CIRCLE ELLIPSE 



TYPE 9 n ELLIPTIC SPIRAL 
WITH REVERSAL 



TYPE 12: LIMIT- CIRCLE 
LIMIT CIRCLE 



TYPE 16: UNBOUNDED LIMIT CIRCLE 
RADIAL AT INFINITY 



TYPE 10: PARABOLIC SPIRAL 



TYPE 12r: LIMIT-CIRCLE 
LIMIT CIRCLE 
WITH REVERSAL 



TYPE 17: UNBOUNDED SPIRAL 
RADIAL AT INFINITY 


Figure 24-Types of orbits that exist in the equatorial plane of the Kerr gravitational field. 
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Table 6— Complete classification of the orbits of a test particle in the equatorial plane of the Kerr gravitational field by region and 

subregion of the master diagrams. 


Schwarzschild Metric 
a = 0 


0<a < 1 


Subregion 3 




1 <a <°° 


Admissible 

Orbits 


It/, Sr, 8, 12 
It/, 5,8, 12 
It/, 5,8, 12r 

It/, 2, 5r, 8, 12 
It/, 5, 9, 12, 13 
lt/,5,9, 12r, 13 
lt/,2,5,8, 12r 
1C, 2, 5, 8, 12 

It/, 4, 5r, 8, 12 
It/, 7, 9, 12, 13 
lt/,7,9, 12r,13 
lt/,4, 5, 12r, 8 
1C, 4, 5, 8, 12 

It/, 7r, 8, 12 
1C, 7, 8, 12 
It/, 7, 8, 12r 
1C, 7r, 8, 12 
1C, 7, 8, 12 
1C, 7, 8, 12r 

It/, 5r, 8, 12 
It/, 5, 8, 12 
It/, 5, 8, 12r 
It/, 5r,8, 12 
It/, 5, 8, 12 
lt/,5,8, 12r 


Subregion 3 


Admissible 

Orbits 


1 U, 5r, 8 la 

lt/,5,8 lb 


It/, 2, 5r, 8 2a 
lt/,5,9, 13 2b 


lt/,2,5,8 

lt/,4,5r,8 3a 
lt/,7,9, 13 3b 


lt/,4, 5,8 


It/, 7r, 8 
1C, 7, 8 


1C, 7r, 8 
1C, 7, 8 


It/, 5r, 8 
lt/,5,8 


It/, 5r, 8 
lt/,5,8 


Subregion 3 


Admissible 

Orbits 

















lid Metric 
0 


0<a< 1 



6a 

6b 

6a* 

6b* 

1 -2 (5 line) 15,9 1 a-2a (5 line) 

lb-2e (5 line) 

1 c-2d ( 5 line) 

2- 5 ((/ line) l£/,5,8 2a-5a ((/line) 

2b-5b ((/line) 
2c-5c ((/ line) 
2d-5c ((/line) 

3- 4 ((/ line) 1 (7, 7 , 8 3a4a ((/ line) 

3b-4b ((/line) 
3c-4c ((/ line) 
3d-4c ((/ line) 

2-3 (IV = 0 line) 3,9 2a-3a (IV = 0 line) 

2b-3b (W = 0 line) 
2c-3c (W = 0 line) 
2d-3d (W = 0 line) 
2e-3e (H 7 = 0 line) 


Admissible 

Orbits 

1(7, 7r, 8, 12 
1(7, 7, 8, 12 
1 (/, 7r, 8, 12 
1 (/, 7, 8, 12 

15, l(/,5r,8, 12 
15,1(7,5,8,12 
15,1(7,5,8, 12r 

1(7,5,8,12, 12r 
1(7,5,8,12 
1(7,5,8,12, 12r 
1(7,5,8,12, 12r 

1(7, 7, 8, 12, 12r 
1(7,7,8,12 
1(7,7,8,12, 12r 
1(7,7,8,12, 12r 

1(7, 3, 5r, 8, 12 
1(7,6,9,12,13 
1(7, 6, 9, 12r, 13 
1(7, 3, 5, 8, 12r 
1(7, 3, 5, 8, 12 


4-5 (W = 0 line) 


10 


4a-5a (W = 0 line) 
4b-5b 0V = 0 line) 
4c-5c (W = 0 line) 


1(7, 6r, 8, 12 
1(7,6,8,12 
1(7, 6, 8, 12r 


ic 


a = 1 

1 <c<°° 

Subregion 4 

Admissible 

Orbits 

Subregion 3 

Admissible 

Orbits 

6a 

1(7, 7r, 8 

6a 

llr 

6b 

1(7, 7, 8 

6b 

11 

6a* 

1(7, 7r, 8 

6a* 

llr 

6b* 

1(7, 7, 8 

6b* 

11 

la-2a (5 line) 

15, 1 (7, 5r, 8 

la-2a (5 line) 

15, 9r 

lb-2e (5 line) 

15,1(7,5,8 

lb-2b (5 line) 

15,9 

2b (S-D line) 

15,9 



2a-5a ((7 line) 

1(7, 5, 8, 12r 

2a-5a ((7 line) 

1(7, 5, 8r 

2b-5b (U line) 

1(7,5,8,12 

2b-5b((71ine) 

1(7, 5, 8 

2b-2e (U-D line) 

1(7, 5, 8 



3a-4a ((7 line) 

1(7, 7, 8, 12r 

3a-4a ((7 line) 

1(7, 7, 8r 

3b4b ((/line) 

1(7,7,8,12 

3b4b ((7 line) 

1(7, 7, 8 

3b-3e ( U-D line) 

1(7, 7, 8 



2a-3a (W = 0 line) 

1(7, 3, 5r, 8 

2a-3a (IV = 0 line) 

3, 9r 

2b-3b (W = 0 line) 

1(7, 6, 9, 13 

2b -3b (IV = 0 line) , 

3,9 

2e-3e (^ = 01ine) 

1(7, 3, 5, 8 



4a-5a (W = 0 line) 

1(7, 6r, 8 

4a-5a (IV = 0 line) 

lOr 

4b-5b (IV = 0 line) 

1(7, 6, 8 

4b-5b (IV = 0 line) 

10 




















4a*-5a* (W = 0 line) lt/,6r,8 
4b*-5b* (W = 0 line) IU, 6, 8 


4a*-5a* (W = 0 line) 
4b*-5b*(JV = 01ine) 


n(2b, 3b, 4b, 5b) 


la-lb (S’ line) 
6b* (E line) 
5a*-5b* (Fline) 


lt/,6,8, 12r n(2a, 3a, 4a, 5a) 
It/, 6, 8, 12 n(2b, 3b, 4b, 5b) 

lt/,6,8 


It/, 8, 14 
It/, 7, 8 
It/, 8, 14 


n(E line, IV = 0 line) 8,16 


2b (H line) 
3b (H line) 


2b (D-S line) 
2b-2e (D-t/line) 
3b-3e (D line) 


It/, 5, 13 
It/, 7, 13 


15,9 
It/, 5, 8 
It/, 7, 8 


la-lb (triine) 

6b* (/; line) 

5a*-5b* (E line) 

n(E line, W = 0 line) 

lb (H line) 

2b (H line) 

3b (//line) 


It/, 6, 8r 
lt/,6,8 


n(2b, 2e, 3b, 3e) lt/,6,8 


a Regions or subregions separated by a hyphen denote the boundary common to the regions or subregions. The symbol n denotes the boundary point common to the 
subregions that follow it. 
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DISCUSSION 


In this section, four topics will be discussed: (1) the determination of direction of orbital motion 
for a given set of constants of motion (h, e ) for the test particle in the Kerr field of a central body with 
mass m = 1 and a given specific angular momentum -a, (2) the scaling of master diagrams and solutions 
from the case in which (3 = 2m = 2 to one in which (3 =£ 2, (3) the conversion of the results obtained in 
the S coordinates to those in the E coordinates, and (4) the behavior of solutions in the neighborhood 
of the pseudosingularity at u = u D . 

Direction of Orbital Motion 


From Appendix D, which illustrates the admissible types of orbits in each subregion of a master 
diagram for the cases of 0 < a < 1 , a = l, and 1 < a < it is observed that the mathematical solution 
u(4>) of Equation 44 has two branches at every point in the 0, u plane in which the one-parameter 
family of solutions F(<f> + c, u) = 0 is defined. One would expect that the branch to be taken is deter- 
mined by the specific angular momentum h of the test particle. Such is the case in the Schwarzschild 
gravitational field. For the present case there is the more general condition given by Equation 219, 


d<j> u 2 E{u) u 2 [~j3(7z + ea)u + h] 

ds D(u) a 2 u 2 - (3w + 1 


( 220 ) 


in which the proper time s is taken to be always increasing. This equation thus determines the direc- 
tion of orbital motion. 


For an illustrative case, consider a central body with (3 = 2 and a = 0.7 and a test particle with 
h < 0 and e < 0 such that (h, e ) lies in region la in the third quadrant of the master diagram given by 
Figure 1 5. (There is also a subregion la in the first quadrant for which h > 0 and e > 0.) Table 6 
indicates that there are four kinds of orbits for this subregion: IU, 5r, 8, and 12, which are then 
obtained from Figure 24. In subregion la of Figure 15, h + ea < 0; and Figure D1 shows that 
0 < itr? < u D1 <u D 2 . With this information, Figure 25 was constructed, which indicates that 0 
decreases for u g3 < u <u E and u Dl < u < u D2 and 0 increases for u £ <u <u D1 and u D2 < u. 
Consequently, the orbits have direction of motion as shown in Figure 26. 


Scaling of Master Diagrams and Solutions 


In this work, the value (3 = 2 was arbitrarily set by taking the central body to have mass m = 1 . 
The results of this study can easily be applied to the case (3^2: Suppose that for a given set (h, e, a, (3) 

there is the solution 

u = w(0; /i, e, a, (3) (221) 

of the differential equation 



D 2 (u) 

= B(u ) . 

E 2 («) 


(44) 
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Next, let a set of new variables (0, u, h, e, a, j?) be defined by 

0 = 0 ' 


_ u 

u~ X 


h = 'kh 


e = e 
a = \a 
0 = X0 

/ 

where X is an arbitrary constant not equal to zero; then, it is obvious from Equations 45 to 47 that 

\ 


( 222 ) 


and that 


B(u;h, e, a, (3) =B(u; h, e, a, j3) 
D(u;a,P)=D(u;a,i 3) 

E(u; h, e, a, = \E(u;h, e, a, j3) 


(223) 


' du\ ? 

l 

idu 

* 

Ft 

X 2 

V0/ 


(224) 


From Equations 44, 223, and 224, it follows that 

u = w(0; h, e, a, (3) (225) 

is automatically a solution of the differential equation 

(226) 

Hence, the transformations in Equation 222 yield the desired scaling of the master diagram and the 

solution. 


(duS* B(u;h, e, a, (3)D 2 (u;a, (3) 
\d4>j E 2 (ii;h,e,a,]3) 


As a concrete example, consider the case of a central body with mass m = 5 (i.e., j3 = 10) and 
specific angular momentum -3 = - 3 and a test particle with constants of motion h = l and 3=11. 
Because the master diagrams in this work are constructed for j8 = 2, from Equation 222 it follows that 



the master diagram for which 



(227) 

(228) 
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and the subregion containing the point 



( 229 ) 


are then used to obtain the various types of admissible orbits and the direction of motion. From these 
orbits in the r, 0 plane, the orbits in the r, 0 plane are obtained from the relations r = 'Kr = 5r and 0 = 0. 
That is, the desired admissible orbits are stretched by a factor of 5. 


Conversion of Results From S to E Form 


The present qualitative solutions are obtained by considering the differential equation that is in the 
S form, Equation 44. However, to obtain the solutions in the E form for which the differential equa- 
tion is 



C* = 0, 


(230) 


where A *, B*, and C* are given by Equations 1 9 to 2 1 , it is not necessary to solve this more compli- 
cated equation. Instead, Equations 10, 12, 14, and 46 are used to obtain 


where D¥= 0. Therefore, 


du* _ du 
d<f>* d<j>* 


du 

< 20 * 


du d<j> 
<20 < 20 * 


a dr*\ 
A d(j>*) 



D du/d<j) 
D + a du/dQ 


(231) 

(232) 


In view of the fact that (h * , e * ) = (h , e) as proved earlier, Equation 23 2 gives the solution u(<p*;h*,e*,a,0) 
in the E form from a knowledge of the solution w(0; h, e, a, j3) in the S form when D i= 0. 

Because <2w/<20 -> ±°° as u -* u E , it is concluded from Equation 232 that 

du D{u £ ) 

Urn = . ' (233) 

duld<t>-+± ~ <20* a 
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Equation 23 1 indicates that if D ¥= 0, then 


du „ du 

— — 0 ■ <==> ' 

<20 < 20 * 


= 0, 


(234) 


which holds at u = u B , as can be easily seen from Equation 44. 

For the case in which <2w/<20 ^ 0 or (2zv/<20 =£ ±°° Equation 231 may also be used to study the 
qualitative nature of the conversion of <2«/<20 to du/d<j>* as follows: For simplicity, let a > 0. Then, 
there are two cases: D > 0 and D < 0. 

For D > 0, there are four subcases. Subcase 1 has the conditions 


^>0 

<20 


so that 


or 


In subcase 2, 


du 

d</>‘ 


> 0 , 


Q< du/d<p* _ a_ du_ < 
du/d<j) D d<i>* 


du du 

d<j>* <20 


du 

(20 


>0 


so that 


dU ^ r, 
— <0 , 
( 20 * 


0> du/dT = ,.£ du_ >0< 

du/d<p D <20* 


which is a contradiction; hence, this case is impossible. For subcase 3, 

^<0 

<20 


du 

c20* 


< 0 , 


so that 


du/d<b* a du 

0 < — = 1 >1 

<2«/<20 D (20* 


(235) 
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or 



du < du 
d(t>* d(j> 

( 236 ) 

Subcase 4 has the conditions 

^<0 

</0 



du ^ „ 

>0, 

<20* 


so that 

du/d<p* a du 

o> - 1 

du/dcj) D <20* 


which is a possibility and for this 

du du 

<20* <20 

( 237 ) 

Similarly, the corresponding four subcases forZ) < 0 are the following: 


(1) 

^>0 



du ^ ^ 

-rz >0; 
<20* 



du du 

<20* <20 

( 238 ) 

(2) 

^>0 

<20 



du „ „ 

<0 ; 

<20* 



<2w <2u 

<20* <20 

( 239 ) 

(3) 

^<0 

<20 



<2« ^ „ 

<0 ; 

<20* 
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( 240 ) 


du du 

d<t>* d<j) 


( 4 ) 


du 

d<f> 


<0 


which is impossible. 


du 

d<{>* 


> 0 , 


Thus, Equations 233 to 240 are sufficient to provide all the pertinent information for conversion 
of any solution in the <p, u plane to one in the 0*, u plane, provided that D # 0 anywhere in the 
solution. Because Equation 230 is also autonomous, the one-parameter family of solutions is given by 


G(<j>* + c, u) = 0 


(241) 


so that all the remaining members in the family can be obtained from any one member by a simple 
translation. Figures 27 and 28 illustrate, respectively, the transformation of «(0) to w(0*) in the 
neighborhood of a turning point and of a reversal point. 



0 * 



u(0) 


u(0*) 


P* 


0 





Figure 27— Conversion of S' form to E form in the neighborhood of a turning point/ 5 . 
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D>0^ 


i D 

/^ slo P e T 




D<0 

> 



Figure 28— Conversion of S form to E form in the neighborhood of reversal point Q. 


Pseudosingularity at u = u D 

In the transformation from E to S coordinates given by Equations 10 to 13, a pseudosingularity 
is introduced when A~0; i.e., D(u) = 0. This means that the transformation of the qualitative results 
back from S to E coordinates holds only if D =£ 0; i.e., u=fru D . 

It remains now to study the behavior of solutions in the neighborhood of u = u 0 . For this, 
we must go directly to the equation of motion in E form (Equation 17) which, in view of u = u*, 
e = e*, and h = h*, may be written simply as 

( du \ 2 du 

Wi) + 2flS(w) ~ B(m)D{u) = 0 

or 

du _ -aB ± yjB{a 2 B + AD) 
d(j>* A 


( 242 ) 

( 243 ) 
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= 0 

U=U D 

2aB(u D ) 
~ a (u d ) 


B(u d ) = 0 






and Equation 244, it follows that only the neighborhood of u = u D =f = u E need be studied because 


otherwise the only possibility is 



( 246 ) 


This equation and the results of the conversion from S to E form suffice to determine the qualitative 
nature of the solution u(0*) in the neighborhood of u = u D = u £ . 

Remembering that by Equation 195 B(u D ) > 0 under these conditions and using Equations 236 
to 241 and 244, the six cases for the solution u(0*) in the neighborhood of u = u D ^u E are deter- 
mined as illustrated in Figure 29. 


Finally, Equation 243 indicates that, at a reversal point in the solution u(tp*), A(u) = 0 must hold. 
In other words, 

(247) 


du 

— - = ±oo => u = u 

d<j>* " 


where 


h 2 +(1 - e 2 )a 2 

u A = , (248) 

|3 (h + ea) 2 


which is obtained from Equations 19, 131, and 134. 


CONCLUSIONS 


The main results are summarized as follows: 

( 1) For equatorial orbital motion, the differential equation in the S coordinates is used because 
of its simple form. 

(2) For j8 = 2 and 0 < a < °°, several master diagrams are obtained for the h, e space of parameters. 

(3) These master diagrams are of three types corresponding to the ranges 0 <a < 1 ,a = 1, and 
1 < a < 

(4) Each master diagram is divided into subregions such that each has the same types of admis- 
sible orbits. 

(5) There are 25 essentially different types of orbits for a =£ 0 (Kerr metric) as compared to 1 1 
for a = 0 (Schwarzschild metric). 

(6) Among the most interesting results are orbits with reversal points or cusps. 

(7) A method is given for scaling master diagrams and solutions for the case of fi =£ 2. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland, January 27, 1972 
630-22-01-01-51 
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Appendix A 


A QUALITATIVE STUDY OF THE COMPLETE SET OF SOLUTIONS OF THE 
DIFFERENTIAL EQUATION OF MOTION OF A TEST PARTICLE IN THE EXTERIOR 
SCHWARZSCHILD GRAVITATIONAL FIELD 


F. K. Chan 

Department of Space Science and Applied Physics 
Catholic University of America 
Washington, D. C. 


INTRODUCTION 

The relativistic orbits of a test particle around a centrally attracting mass have been studied time 
and again since Einstein (References A1 and A2) first studied the quasi-elliptic orbits in connection 
with his analysis of the precession of the orbit of Mercury around the Sun. Later work on this type of 
orbit has been contributed by several authors among whom are Droste (References A3 and A4), 
de Sitter (References A5 and A6), Greenhill (Reference A7), and Forsyth (Reference A8). Although 
Droste and Morton (Reference A9) have also considered various other types of orbits, their works are 
not really complete because they do not give an exhaustive account of all possible types of orbits. It 
was not until the work of Hagihara (Reference A 10) that a truly complete and exhaustive analytical 
treatment of the whole problem was first published. Work following this includes that of Kaplan 
(Reference All) who studied only the circular orbits, Darwin (References A1 2 and A 13) who obtained 
a fairly detailed but not totally complete account of the various types of orbits, Metzner (Reference 
A 14) who essentially completed the studies of Darwin, and Galkin (Reference A 15) who made some 
study, though not elaborately revealed, on certain aspects of relativistic orbits. A fairly complete 
account is also given by Arzelies (Reference A 16). However, the works of Hagihara, Darwin, and 
Arzelies, to mention the most authoritative sources on the subject, are rather involved and rely heavily 
on the use of elliptic functions. A different and apparently much simpler approach to the analysis of 
these orbits, involving a qualitative study of the solutions of the governing differential equation of 
motion of the test particle, was obtained by the present author in 1965 with only a knowledge, at that 
time, of the first paper of Darwin. The analysis given here is essentially taken from his unpublished 
work on certain aspects of nonlinear ordinary differential equations in which he studied, as an exam- 
ple, the relativistic differential equation of motion of a test particle in an exterior Schwarzschild 
gravitational field. 
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ANALYSIS 


The motion of a test particle in the exterior Schwarzschild gravitational field produced by a 
spherical body of mass m is given by the well-known equation 

/dll's} 2m W 

[ — ] = 2 mu 3 - u 2 + — u + — , (Al) 

\d<j)/ h 2 h 2 

where h and W are constants of the motion dependent only on the initial conditions of the test parti- 
cle. The quantity h may be interpreted as the angular momentum per unit rest mass of the particle, 
except now, as contrasted with the Newtonian analog, it is defined in terms of the rate of change of 
the angle <f> with respect to the proper time s of the particle, i.e.. 



(A2) 


The quantity W may be interpreted as a measure of the “energy” of the particle and in the Newtonian 
limit it reduces to 2 E, where E is the total energy per unit mass of the particle. 

Topological properties of the solutions satisfying Equation Al will be studied by considering the 
geometric-topological nature of the surface u, v ) = 0, where for convenience ( du/dq i>) 2 is denoted 

by v. In Equation Al, <p does not appear explicitly; consequently, in the <p, u, v space, the surface 
<1>(<P, u, u) = 0 is a cylinder with a cubic curve as the trace in the u, v plane. Because the radius r is 
defined to be nonnegative, the region to be studied is in the <t>, u plane for which u > 0. Furthermore, 
it is necessary that v>0; therefore, the region to be considered in the 0, u plane is that for which 
u > 0 and v > 0. Consequently, conditions are sought for which there are zero, one, two, or three 
positive roots of the trace cubic curve. A general approach in this direction would require the algebraic 
knowledge of roots of the cubic equation and this is available; but the analysis becomes very cumber- 
some algebraically. In addition, it becomes more involved, for instance, in the case of one real posi- 
tive root whether the remaining two roots are real but negative or are complex conjugates. 

However, in the following analysis the surface <k(<£, u, v) = 0 will be studied without having to 
consider these difficulties. In this process, a geometric meaning will be given to the angular momen- 
tum k and the energy W. Also, because the surface is cylindrical, only the trace in the u, v plane will 

be considered. 

For convenience, the expression on the right-hand side of Equation Al is denoted by /(«), so that 
(duS 2 2m W 

I — 1 = v 2 = f(u) = 2 mu 3 - u 2 + u + — , (A3) 

\d(j>/ h 2 h 2 

2m 

f (u) = Emu 2 -2 u+ — , (A4) 

h 2 

f"{u) = 12 mu - 2 . (A5) 


78 



It is seen that f”(u) does not involve h or IV and is therefore independent of the initial conditions of 
the test particle. In other words, the curve f"(u) is the same for all orbits in the field of a gravitational 
mass m; it is a straight line intersecting the u axis at u = 1 /6m or r = 6m. 


The first integral of f”(u), i.e., /'(«), is a parabola with the axis of symmetry at u s = l/6m. From 
Equation A4, it is seen that the angular momentum appears in the constant of integration. If they 
exist, the intercepts u'^ and u'^ 2 of this parabola f\u) with the u axis are given by 


A,f2 



(A6) 


which indicates that 0 < u'^ < l /6m and 1 /6m < u^ 2 < 1/3 m. A sketch of this family of parabolas, 
with h 2 as the parameter, is shown in Figure Al. 

The second integral of f"(u), i.e.,/(u), is a cubic curve and involves a second arbitrary constant 
given by W/h 2 . That is, the quantity W which is a measure of the energy of the particle appears only 
in the second integral of/"(w), whereas the angular momentum h makes its first appearance in the 
first integral of f'(u). A sketch of this family of cubic curves, with W as the parameter, corresponding 
to the four cases of the parabolas given in Figure Al, is shown in Figure A2. 


The two intercepts, if they exist, of the parabola /'(«) determine the maximum and the minimum 
of the cubic curve /(u) and are in turn determined by the angular momentum. Consequently, the role 
of the angular momentum h is to determine the distance of separation between the points of maximum 

and minimum in the /(«) curve, whereas the 
quantity W or rather W/h 2 determines its verti- 
cal displacement. 



curves. 


The region in the 4>, u plane for which u > 0 
and v > 0 is a simply connected one for the case 
in which 0 < h < 1 2 m 2 , as seen in Figure A2. 
However, for h 2 > 12m 2 , it is possible for this 
region to consist of two simply connected re- 
gions that are disjoint from each other (see Fig- 
ures A2c and A2d); this situation corresponds to 
the case in which the cubic curve f(u ) has two 
or three real positive roots. 

Attention will be focused on the trace cubic 
curve as shown in Figure A2c because it con- 
tains features that are shared by the other three 
types. For example, consider the surface 
<i>(0, u,v) = 0 with trace on the u, v plane as 
shown in Figure A3. The quantities u'^ and u'^ 2 
(«j and Wjj in the figure) correspond to the two 
roots of the parabola /'(«) = 0 and are given by 
Equation A6, whereas the quantities , iy 2 , 
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v or w 



Figure A3— Trace cubic curve of class c with 1/1/ < 0 and 1/4/77 < u' u < 1/3 m. 


and Uf 3 (Wj, « n , w m ) correspond to the three roots of the cubic curve /(u) = 0. From this, the surface 
F(0, w, w) = 0 is obtained (Figure A3) where w = du/dctr, this surface corresponds to Equation A1 in 
the 0, u, w space. 


From F(0, u, w) = 0, it is easily seen that in the 0, u plane the three lines u - , u^ 2 , and u^ 3 

correspond to envelopes or singular solutions; this result follows from the sufficient conditions (see 
Reference 17): 

F(0, u,w) = 0 , (A7) 



9 F 9 F 

— + w — = o , 

90 bu 


(A8) 

(A9) 
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From the trace of the surface F{<j>, u, v) = 0 as shown in Figure A3, the family of solutions u = u(<p) as 

shown in Figure A4 is obtained. 

The solutions u < u(<p) < u f2 are periodic and almost sinusoidal. In fact, they correspond to 
the elliptic orbits of the Newtonian theory for which the trace in the u, v plane is a parabola so that 
the corresponding trace in the u, w plane is a circle. As contrasted with the Newtonian theory, the 
relativistic theory gives an additional region u(<p) > w^ 3 , where solutions exist. 

Let us now study the effects of shifting vertically the surface <3?(0, u, v) = 0 or, equivalently, the 
trace in the u, v plane. This corresponds to the parameter W assuming various values. 

If this cubic curve is lowered, the u{4>) solutions as shown in Figure A4 retain their topological 
properties until the situation is reached where the roots and u^ 2 coincide, given a double root 

u = Uj- l = u^ 2 = u'jy . Here the whole family of solutions < u(0) < u^ 2 has degenerated into a single 

line u = u'f X which is now an “isolated solution” and corresponds to the case of circular orbits. (Per- 


U 



Figure A4— Solutions u(0) corresponding to differential equa- 
tion given by trace cubic curve of Figure A3. 


haps it should be pointed out here that the three singular solutions u = Up, Up , and u y 3 as shown in 
Figure A4 are not physically admissible as real orbits because they are three circles of different radii 
but with the same angular momentum h and also the same “energy” measure W. However, they are 
mathematically admissible because they are also solutions of the differential equation of motion.) On 
shifting the trace curve lower, we see that only the region u(<p) > Up continues to persist and nothing 
essentially different happens thereafter. 

Now if the trace curve as shown in Figure A3 is raised, again there are no topological changes in 
the solutions u(<p) until either one of two situations occurs: Either the cubic curve at u - 0 becomes 
positive before the minimum point at u = u'p does or vice versa; i.e., either /(0) > f(u'p ) or vice versa. 
As shown later in this analysis, the criteria for these two cases are l/4m < u'p < l/3m and l/6m < 
u 'f 2 ^ respectively. However, for the present discussion, it does not matter which happens first, 
and only the former case will be considered in detail. When this situation is reached, there is no change 
in the mathematical properties of the solutions u(0); they are still oscillatory in nature except that now 
they assume at various points the value u((p) = 0. Physically, however, it means that the particle start- 
ing at some initial radius greater than 1 /« either begins to move inward, corresponding to one family 
of solutions, to some smaller radius and then flies off to infinity or, on the other hand, it begins to 
move outward to infinity, corresponding to the other family of solutions. In either case, the particle 
“barely” makes it to infinity. 

Now, for the sake of convenience, an orbit is classified as elliptic-, parabolic-, or hyperbolic-like 
according to whether W is less than, equal to, or greater than zero. This criterion reduces to that for 
Newtonian orbits except that now the orbits are not true ellipses, parabolas, or hyperbolas. Also, 
the orbits shown in the upper region of Figure A4 will be referred to as “captured” orbits because they 
are bounded and eventually fall into the central attracting mass because u($) -* °° as 4> -*■ ± °o 


If we continue to raise the trace cubic curve higher, no essential changes occur in the orbits ex- 
cept that they become hyperbolic-like. Eventually the situation is reached in which the roots Up and 
Up coincide, thus giving a double root u = Up = Up - u'p. The traces on the u, v plane and the u, w 
plane now appear as shown in Figure A5. 


If the variable r] is defined as 


n = u- u'p , 

then it may be verified that Equation A1 can also be written as 

/ du \ 2 , , 

(All) 

U/ = " 277 +flsT? ’ 

where 

(A 12) 

a 2 = 1 2 mu'p - 2 

(A 13) 

fl 3 = 12m . 

(A 14) 

Because h 2 > 12 m 2 for the type of trace cubic curve being considered (see Figure A2c), from Equa- 
tions A6 and A13 it follows that we have a 7 > 0. 
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Figure A5— Trace cubic curve of class c with l/V> 0 and 
1/4 m<u' u < 1/3 m. 


In a small neighborhood of the line u = u^ 2 in the 0, u plane 

du \ 2 , 

d<t>) 2 


or 


d V\ 2 


(A 15a) 
(A 15b) 


from which the whole family of solutions is given by 

r\ = C exp \/a ~2 0) (A 16a) 

or 

u = u'f 2 + C exp (±y/a^cp^ . (A16b) 

In other words, instead of having the regular type of envelope at u = u '^ 2 , there is an “asymptotic” 
envelope and the family of solutions w(0) is shown as sketched in Figure A6. For completeness of 
discussion, if k 2 = 12m 2 , thena 2 = 0 and Equation A12 becomes 

2 

= a 3 v 3 , (A17) 

from which is obtained 

u = u'f 2 + 4 (C ± sj~a 3 2 , (A 1 8) 

thus giving solutions that are essentially the same in nature as those sketched in Figure A6. 
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Figure A6— Solutions u((j>) corresponding to differential equation 
given by trace cubic curve of Figure A5. 


Corresponding to the two solutions passing through a given point in the 0, u plane, from Figure 
A 6 is seen that a test particle, starting at some initial radius greater than l/u| 2 , would either spiral in 
toward a circle of radius 1 /u'^ 2 or would fly off to infinity. Similarly, if it is initially inside that circle 
it would either spiral out toward it or would fall into the center of the attracting mass m. The circle 
of radius 1 /w| 2 m this case is a physically admissible orbit and any particle starting to move on it will 
continue to do so. However, this circular orbit is unstable with regard to any small perturbation in 
the sense that any outward disturbance would send the particle off to infinity whereas any inward one 
would send it into the center of the attracting mass. In contrast, the circular orbits previously dis- 
cussed with radius \/u ^ are stable because any perturbation, inward or outward, would result in an 
elliptic-like orbit close to the unperturbed circular one. 

Considering now the case in which the minimum point of the trace cubic curve becomes positive 
before the point at u = 0 does (i.e.,/(w| 2 ) > /( 0)), it is seen that there are essentially no differences 
from the previous case as illustrated in Figures A5 and A6 except that now the orbits are bounded 
instead of unbounded, as shown in Figures A7 and A8. 

Summarizing the more important results obtained up to this point, the conditions 0 < u ^ < 1 /6m 
and l/6m < < l/3m have been established and, from the nature of the solutions illustrated in 

Figures A4, A6, and A8, it is seen that circular orbits with radius 6m < r < °° are always stable whereas 
those with radius 3m < r < 6m are always unstable with regard to perturbations. Moreover, for 
4m < r < 6m the outward perturbed orbits are bounded whereas for 3m < r < 4m they are unbounded; 


85 



in both cases the inward perturbed orbits eventually fall into the central attracting mass m. The main 
result is that no circular orbit can exist with a radius r<3m. It has also been established that, be- 
cause l/6m < u'f 2 < 1/4 m corresponds to the case of/(0) </(w| 2 ), all the elliptic-like orbits must 
have a perihelion not less than 4m. (See Figure A7.) Similarly, because 1/4 m < w/ 2 < l/3m corre- 
sponds to the case of/(0) >/(w| 2 ), therefore all hyperbolic-like orbits must have a perihelion not less 
than 3m. (See Figure A5.) Also, the minimum perihelion for a parabolic-like orbit is 4m. 

If both the trace cubic curves as shown in Figures A5 and A7 are raised, then the two correspond- 
ing families of solutions u(<p) are obtained as shown in Figures A9 and A 10. On further raising the 
trace cubic curve in Figure A7 so that it is everywhere positive for u > 0, then the family of solutions 
in Figure A10 will look essentially like that in Figure A9. No further changes occur in either of these 
two cases on continued vertical displacement of this trace curve. 

The orbits in the lower region of Figure A6 (W > 0) are referred to as “hyperbolic limit circles” 
and the corresponding ones in Figure A8 ( W < 0) as “elliptic limit circles.” The orbits in the upper 
region in these two diagrams are called “limit-circle spirals.” For the case of W = 0 in Figure A6 or 
A8, orbits occur in the lower region which are referred to as “parabolic limit circles.” Finally, those in 
Figures A9 and A10 are called “hyperbolic spirals” and “elliptic spirals,” respectively. As contrasted 
with the mathematical representation of these various types of orbits in the <j>, u plane, Figure A! 1 
illustrates the physical appearance of these orbits. In this diagram, the name given to describe each 
orbit attempts to illustrate its properties at r v and r. . Thus, elliptic spirals and captured orbits 
are classified as of the same type. 



Figure A9— Solutions u{<j>) corresponding to differential 
equation given by trace cubic curve of Figure A5with a 
larger positive value of W. 


Figure A10— Solutions u(4>) corresponding to differential 
equation given by trace cubic curve of Figure A7 with a 
less-negative value of W. 


87 



TYPE Is: STABLE CIRCULAR ORBIT / ^ ^ 

TYPE lu: UNSTABLE CIRCULAR ORBIT TYPE 2: ELLIPTIC-LIKE ORBIT TYPE 3: PARABOLIC-LIKE ORBIT TYPE 4: HYPERBOLIC -LIKE ORBIT 

(See Figures A13 and A14) (See Figure A4, W < 0) (See Figure A4, W = 0) (See Figure A4, W > 0) 



TYPE 5: ELLIPTIC LIMIT CIRCLE TYPE 6: PARABOLIC LIMIT CIRCLE TYPE 7: HYPERBOLIC LIMIT CIRCLE TYPE 8: LIMIT-CIRCLE SPIRAL 

(See Figure AS, W < 0) (See Figure A6, W = 0) (See Figure A6, W > 0) (See Figures A6 and A8) 


C3 

TYPE 9: ELLIPTIC SPIRAL 

(See Figures A10, W < 0) (See Figure A9, W = 0) (See Figure A9, W > 0) 

Figure A1 1— Physical appearance of various types of orbits. 

Although the trace curve in the u, v plane as shown in Figure A2c has been discussed specifically, 
the cases as shown in Figures A2b and A2d in the limits of u'j X = u/ 2 - l /6m and w^ 2 = l/3m, respec- 
tively, have actually also been covered. The remaining case as shown in Figure A2a is virtually the 
same as that of Figure A2b as far as the family of solutions w(0) in the (j>, u plane is concerned. 

The next step is to show that u’ f2 = 1 /4m is the line of demarcation for the two cases described 
earlier; i.e., whether/(0) > /(wT ) or vice versa. To do this, consider the following theorem: 

For a general cubic curve of the form 

y = ax 3 + bx 2 + cx + d , 

let x = x be the axis of symmetry of the cubic curve, let x Q be the value of x corresponding to either 
the local maximum or the minimum value of y, and let x Q be the value of x such that 

y(x Q ) = y(x Q ) 

forx 0 Then the following relationship is always true: 

2(x 0 - x) = x - x 0 . 

Accepting this theorem and returning to the problem, it can be seen that because the trace cubic 
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curve always has u = I /6m as its axis of symmetry, the following relations result: if 

, = J_ 

f 2 Am 


1 , 1 

< u' < — 

Am f 2 


5m 


then, respectively, 


as shown in Figure A 12. 


J_ , _L 

6m f 2 Am 


/(0) = f(u f2 ) , 

m >f(u’ f2 ) , 
/CO) <f(u' f2 ) , 


The relationship between the radius and the angular momentum (per unit mass) for circular orbits 
is determined next. From the preceding discussion it is noted that circular orbits, whether stable or 
unstable, correspond to the situation in which there is a double root in the trace cubic curve in the 
u, v plane; i.e.,/(wy) = 0 and f'(u^) = 0. For a specific value of angular momentum, these double 

roots are given by Equation A6. Denoting 
the inverse radius of a circular orbit by u c , 
the following relationship is obtained between 
the radius (r = \/u c ) of a given circular orbit 
around a gravitating mass rn and the angular 
momentum per unit mass h which the parti- 
cle must have to remain at that radius: 



1 


6m 



(A 19) 


Figure A 12— The genera! cubic curve illustrating the 
relation x - x 0 = 2{x Q - x). 


It follows that for a given angular momentum, 
there exist two possible circular orbits; in 
Equation A 19, the plus sign is associated with 
unstable circular orbits and the minus sign 
with stable circular orbits. In either case, it 
follows from this equation that circular orbits 
can only exist if the following inequality is 
satisfied: 

12m 2 

V= l >0 . 

h 2 
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In Equation A19 the larger root corresponds to the family of unstable orbits as shown in Figure A13 
whereas the small root corresponds to the family of stable orbits as shown in Figure A 14. 


Next, the relationship between the radius and the “energy” W c for circular orbits is obtained 

from Equation A3: 


/ 2m 

W c =~h 2 \ 2mul - u\ + ~ u c 


(A21) 


Substitution of u c from Equation A 1 9 yields 

1 h2 

S=W. = (1 - a) + (a 3 - 3a + 2) 

3 108m 2 

for 0 < u. < 1 /6m, and 

1 h2 

U=W ,, = ( 1 + ck) (ce 3 - 3a- 2) 

3 1 08m 2 


for 1 /6m < w < l/3m, where 


a= y/V = 



(A22) 


(A23) 


(A24) 


From Equations A22 and A23 it is observed that the continuity result S = U holds at u c - 1 /6m (i.e., 

a = 0) as expected. 


f(u) 



Figure A1 3— Criterion for unstable circular orbits, 
3m 6 m. 


U 


Figure A14— Criterion for stable circular orbits, 
6/7? < r < °°. 
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Introducing, for convenience, the variable £ defined by 

h 2 

|= , (A25) 

m 2 


the following asymptotic representations for S and U when £ -> °° are obtained: 



(A26) 


1 2 1 

£ 7 = — £ +0 - 

27 5 3 £ 


(A27) 


The complete quahtative nature of the solutions of Equation A 1 for all physically admissible 
ranges of the two parameters h and W has been studied. It is to be noted that the physically admissible 
ranges of h and W are 0 < h 2 <°° and - 1 < W < the latter condition resulting from the following 
equations: 


1 - 


2m\ dt 
ds 


C 


(A28) 


W = e 2 - 1 . (A29) 

All the results are represented in the master diagram shown in Figure A 1 5, which gives the topological 
genera of the orbits corresponding to a given set of values of these two parameters. This diagram is to 
be used in conjunction with Table A1 and Figure A1 1 for determining the physical appearance of 
admissible orbits corresponding to a given set (h, W) in Equation Al. 

Characteristics of the cubic curve f{u) in various regions of the master diagram (Figure A 15) are 
illustrated in Figure A 16 and tabulated in Table A2. 

From the Newtonian theory, the following relationships between the radius (r = l/u c ) and the 
angular momentum per unit mass h and the energy per unit mass E hold for circular orbits: 

h 2 = (A30) 


G 2 m 2 

2E = , (A31) 

h 2 

so that the corresponding, but much less complicated, master diagram appears as in Figure Al 7. Be- 
cause W -*■ 2E in the Newtonian limit, the quantity 2E is used as the variable on the vertical axis. Also, 
the gravitational constant G, which is set to unity in the relativistic diagram, now makes its appear- 
ance in the variable on the horizontal axis. 
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Figure A15— The master diagram for relativistic orbits in terms of W and h, both per unit rest mass of test 

particle. 
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Table A 1— Complete classification of the orbits of 
a test particle in the equatorial plane of the 
Schwarzschild gravitational field. 


Region of Master Diagram 

Admissible Orbit Nos. j 

1 

9 1 

2 

2,9 j 

3 

4,9 

4 

It 

4 * 

il 

5 

9 

5* 

9 

1-2 (S line) 

15,9 

2-5 (U line) 

W,5, 8 

3-4 (U line) 

li/,7,8 

2-3 (W ~~ 0 line) 

3,9 

4-5 (W = 0 line) 

10 

4*-5* (V/ = G line) 

1G 

Intersection of 2, 3, 4, 5 

6,8 

i 



Figure A 16— Characteristics of the f{u) curve in various regions. 


/ 


U 

: W >0 ,V>0 


/ 

y 


u 

; W < 0 , V > 0 
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Table A2— Characteristics of the f{u ) curve for 
various regions of the master diagram. 


Region 


/( w /2> 

D 

w 

u s 

1 

<0 

<0 

>0 

<0 

>0 

2 

>0 

<0 

>0 

<0 

>0 

3 

>0 

<0 

>0 

>0 

>0 

4 

>0 

>0 

>0 

>0 

>0 

4* 

( a ) 

( a ) 

<0 

>0 

X) 

5 

X) 

>0 

>0 

<0 

X) 

5 * 

( a ) 

( a ) 

<0 

<0 

X) 


^u'f \ and u'^2 do not exist for V < 0. 



Figure Ml — The master diagram for Newtonian orbits in 
terms of E and h, both per unit mass of test particle. 


CONCLUSIONS 

From a comparison of the results obtained from these two theories, it is concluded that— 

(1) The S curve is identical with the 2 '£ curve for large values of the variable h 2 /m 2 . 

(2) There is no lower limit to the quantity 2 E for the Newtonian orbits; but one does exist for 
for the energy measure W for the case of relativistic orbits. 













(3) For a given set of {h, 2 E) there exists only one type of Newtonian orbit; for a given set of 
(h, W) there can exist one, two, or even three types of relativistic orbits. 

(4) It is well known that the shape of Newtonian orbits is completely specified by only one 
parameter, the eccentricity e, which is defined by the equation 



However, in contrast, relativistic orbits of a given type cannot be specified by just one 
parameter. 

(5) The relativistic orbits experience no singular behavior at r = 2m. This is an indirect proof of 
the lack of any intrinsic geometrical singularity at the Schwarzschild radius. 

(6) In Newtonian theory, circular orbits are possible for all radii. However, in relativity theory, 
no circular orbit can exist with a radius r < 3m. Moreover, circular orbits are stable if 
6m<r < °° and unstable if 3m < r < 6m. Finally, if 3m < r Am, the outward perturbed 
orbits are unbounded; whereas if Am < r < 6m, they are bounded. In both cases, the in- 
ward perturbed orbits eventually fall into the central attracting mass. 

(7) In relativity theory, elliptic-like and parabolic-like orbits must have a perihelion not less 
than Am whereas hyperbolic-like orbits must have a perihelion not less than 3m. However, 
no such limits exist for their Newtonian analogs. 
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Appendix B 


DETAILED DERIVATION OF RELATIONSHIPS USED IN 
ANALYSIS OF EQUATORIAL ORBITS 


*RQGF THAT B(u' 8 ) > 0 IN REGION 6 


The expression for B(u) is given by Equation 45, which is 

B(u ) = 0(/z + ea ) 2 u 3 - [k 2 + ( 1 - e 2 }a 2 ] u 2 + j3u + e 2 - 1 , 
and the roots of B’(u) - Q are given by Equation 142: 

U Bl,B2 =U s( l ± V^) > 
h 2 +(1 - e 2 )a 2 


where 


“» = 


30(Ji + ea} 2 


V= I 


30 2 {h + ea) 2 
lh z +(1 - e 2 )a 2 ] 2 


1 


= 1 


3 u 2 {h + ea) 2 


In region 6, 


* s <0 : 


0 < F< 1 , 

e 2 - 1 > 0 . 

Substitution of Equations 142, Bl, and B2 into Equation 45 leads to 


B(u' 3 ) = 2p(l±^/V)u s 


- - ( h+ea) 2 u 2 


1 •> 

+ -• §u s + ( e 2 - 1 ) . 


(45) 

(142) 


(El) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 


From Equations B2 and B4 it follows that 



Hence, from Equations B3, B4, and B7, it is seen that 


2/3 (l ± yf v)u s - - (h + ed) 2 u 2 > 0 . (B8) 

Next, consider the remaining group of terms in Equation B6. For convenience, the following defini- 
tion is made: 

0 u 

Fie, h)= — + e 2 - 1 . (B9) 

3 

If Equation B1 is substituted into this equation, then 

F{e, h ) = ih 2 [9(e 2 - 1) + 1] + 18 ae{e 2 - \)h + 9a 2 e 2 {e 2 - 1) + (1 - e 2 )a 2 \ , (BIO) 

9(h + ea) 2 

which is quadratic in h whose discriminant d is given by 

d= [18 ae(e 2 - l)] 2 - 4[9(e 2 - 1) + 1 ][9a 2 e 2 (e 2 - 1) + (1 - e 2 )a 2 ] 


= -3 2a 2 {e 2 - 1) . 

(Bll) 

Because e 2 - 1 > 0 by Equation B5, then 


d< 0 

(B12) 

and so 


B(u' b )*0. 

(B 1 3) 

Next it is noted that 


1 - e 2 


Fie, 0) = e 2 - 1 + > 0 , 

(B14) 


9e 2 

so that, in view of Equation B 1 2, it follows that 

Fie, h) = e 2 - 1 + >0. (B15) 

3 

Hence, from Equations B8 and B15, 

B(u' b )> 0 (B16) 

in region 6. 


PROOF TIHAT/1(u e ) = 0 AND D(u £ ) = 0 SIMULTANEOUSLY ARE IMPOSSIBLE 

From the equation 

D(u ) = a 2 w 2 - J3w + 1 

and the identity E(u £ ) = 0 or, equivalently, 

- h 
E 1 3(h + ea) 


(46) 

(206) 
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the following is obtained: 


D{u e ) = a z u E - fiu E + 1 

a 2 h 2 0h 

= t 1 

fi (h + eaj 2 (3(h + ea) 


The equation 
and Equation 206 yield 

Given the equations 


1 

= ( a 2 h 2 + afieh + a 2 fie 2 ) . 

fi(h + ea) 2 

A(u) = h 2 + (1 - e 2 )a 2 - ft (h + ea) 2 u 

A(u e ) = h 2 + (1 - e 2 )a 2 - h(h + ea) 

= a 2 - a 2 e 2 - hea . 


A(u e ) = a 2 - a 2 e 2 - hea 
= 0 

1 

D(u e )= (a 2 h 2 + afi eh + a 2 fi e 2 ) 

fi (h + ea) 2 


= 0 , 

then substitution of hea from Equation B19 into Equation B20 leads to 


1 

( a 2 h 2 + fia 2 ) = 0 , 

fi ( h + ea) 2 

which is impossible because )3 =¥= 0 and a ¥= 0. 


PROOF THAT A'(u £ ) = 0 IS IMPOSSIBLE 
From the equation 

A(u) = h 2 + ( 1 - e 2 )a 2 - ${h + ea) 2 u , 

it follows that 

A'{u) = -&(h + ea ) 2 

for all u. From this equation and 

h 

UjE (3(h + ea) 

it is seen that 


for any finite u £ . 


A'{u e ) =£ 0 


(B17) 

(48) 


(B18) 

(B19) 


(B20) 


(B21) 


(48) 

(B22) 

(206) 

(B23) 
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SOLUTION OF D(u d } = 0 5 D’(u d } = 0, E{u D ) = 0, AND B"(U D ) = 0 

From Equations 45 to 47 and 190, the following equations are obtained: 

D{u d ) = a 2 u\ - j 3 u d + 1=0, (193) 

D'(u D ) = 2a 2 u D -0 = 0, 


E(u 0 ) = -0(/? + ed)u D + h - 0 , 


B"(u D ) = 60(/z + ea) 2 u D - 2{k 2 + (1 - e 2 )a 2 } = 0 . 
Equations 193 and B24 indicate that u D is a double root of D{u) = 0, so that 

1 

u l = ~ 


or 


Equations B24 and 19 Id yield 


U = + - 

U D ~ 


0 = ±2a 


F 19 Id) 
(191 c) 


= 2 . 


Hence, substitution of u D and 0 from Equations 1 9 1 d and 191c into Equation B25 leads to 


h = - 2ea . 


Then, from Equations B26 to E28, the following are obtained: 

3 pe 2 a 2 u D - (3 e 2 a 2 + a 2 ) = 0 , 

Se 2 a 2 - he 1 a 2 - a 2 = 0 , 


Consequently, Equation B28 then gives 


3e 2 - 1 = 0, 


e 




(191a) 

(191b) 


FACTORIZATION OF B{u £ ) 


From Equations 178, B17, and B18, the following definitions result: 

1 

B(u £ ) = - — A(u £ )D{u e ) , 
cr 
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(B24) 

(B25) 

(B26) 


(B27) 

(B28) 


(B29) 



(B30) 


6 2 eh 


D{u e ) = — 

j3 2 (h + ea) 2 

A{u e ) = -a(eh + ae 2 - a) . 
Equation B30 can be written in the form 


+ P 2 e 2 l , 


D{u~) = 


h 2 


From the definition 


p 2 (h + ea ) 2 

a 2 Up - 1 3u D + 1 = 0, 

the sum (u D x +u D2 ) an( * product u D1 u D2 of the roots are 

J3 

U Dl + U D2 ^ 
a~ 


+ haPe (—\ + a 2 p 2 e 2 ! — j] . 


V -2/ 

\u / 
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u Di U D2 


From this, Equation B32 can be written as 


D(u f ) 


P 2 (h + ea) 2 
a 2 

ft 2 (h + ea) 2 
a* 


[h 2 + hafk(u Dl +u D2 ) + a 2 0 z e 2 u nx u D2 ] 


did: 


(h + a$eu D , ){h + afieu D2 ) 
| h + /H ! 


1 3 2 ( h + ea) 2 \ au 


D 2 / 


au 


DIJ 


Hence, it follows from Equations B29, B31, and B32 that B(u E ) can be written as 

/ 
a 


3(u e ) = — 

(j 2 (h + ea ) 2 

Similarly, from Equations B29 to B31, 


i e & \ ( e ‘ 6 

{eh - a + ae 2 ) ( h H 1 { h + 


au 


Dl/ 


au 


D 2 / 


ae 


B(u f ) = 


p 2 (h + e) 2 


(eh - a + ae 2 ) 


(h\ 2 p 2 h 1 

-) + — - + P 2 

\eJ a e 


J' 


( 200 ) 
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(B34) 


(B35) 


101 




Appendix C 

MASTER DIAGRAMS FOR 0 = 2 AND VALUES OF a RANGING FROM a = 
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Figure C4~Master diagram for a = 0.3, 0 = 2.0 




















6a 


a 



ER UNIT MASS 
for a = 0.5, (3 = 2.0. 

























H U 



114 







SSUW ilNfl U3d ADU3N3 


115 


Figure Cl 1-Master diagram fora = 1.0, fi = 2.0. 
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Figure C13— Master diagram for a = 1.2, 0 = 2.0. 
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Figure C14— Master diagram for a = 1,3, 0 = 2.0, 
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Figure CIS-Master diagram for a = 1.4, 0 = 2.0. 
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ANGULAR MOMENTUM PER UNIT MASS 












SSHW llNfl U3d JL0U3N3 


130 


ANGULAR MOMENTUM PER UNIT MASS 
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ANGULAR MOMENTUM PER UNIT MASS 
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Figure C37— Master diagram for a = 90.0, (5 = 2.0, 





Appendix D 


QUALITATIVE SOLUTIONS IN REGIONS OF MASTER DIAGRAMS FOR 

0<3<1,a = 1, AND 1 <3 <°o 
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Figure DIO— Qualitative solutions on D * line for 0 <a < 1. 
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Figure DIO (concluded)— Qualitative solutions on D 1 line for 0<a < 1. 
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i— Qualitative solutions on E line for , 
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Figure D30-Qualitative solutions in region 1 for 1 <a <°°. 


180 











SOLUTION u {*) 


SOLUTION u (*) 



REGION 4a-PHYSICAL ORBITS 



REGION 4b-PHYSICAL ORBITS 


Figure D33— Qualitative solutions in region 4 for 1 <3 < 00 
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Figure D35— Qualitative solutions in region 5 for 1 < a < 
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Figure D36— Qualitative solutions in region 5* for 1 <a<°°. 
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Figure D37— Qualitative solutions in region 6 for 1 <a <°°. 
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Figure D42 (concluded)— Qualitative solutions on 
U line for 1 <a <°°. 
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